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PREDICTION  OF  VORTEX  SHEDDING  FROM 
CIRCULAR  AND  NONCIRCULAR  BODIES  IN 
SUBSONIC  FLOW 

Michael  R.  Mendenhall 
Daniel  J.  Lesieutre 

Nielsen  Engineering  &  Research,  Inc. 

SUMMARY 

An  engineering  prediction  method  and  associated  computer 
code  VTXCLD  to  predict  nose  vortex  shedding  from  circular  and 
noncircular  bodies  in  subsonic  flow  at  angles  of  attack  and  roll 
are  presented.  The  axisymmetric  body  is  represented  by  point 
sources  and  doublets,  and  noncircular  cross  sections  are  trans¬ 
formed  to  a  circle  by  either  analytical  or  numerical  conformal 
transformations.  The  lee  side  vortex  wake  is  modeled  by  discrete 
vortices  in  crossflow  planes  along  the  body;  thus  the  three- 
dimensional  steady  flow  problem  is  reduced  to  a  two-dimensional, 
unsteady,  separated  flow  problem  for  solution.  Comparison  of 
measured  and  predicted  surface  pressure  distributions,  flow  field 
surveys,  and  aerodynamic  characteristics  are  presented  for  bodies 
with  circular  and  noncircular  cross  sectional  shapes. 

INTRODUCTION 

Current  missile  applications  reguiring  high  aerodynamic 
performance  can  involve  noncircular  body  shapes  in  subsonic  flow 
at  high  angles  of  attack  and  nonzero  roll  angles.  The  angle  of 
attack  range  may  be  sufficiently  high  to  cause  formation  of  body 
separation  vortices,  and  the  body  vortex  shedding  characteristics 
are  directly  influenced  by  the  body  cross-sectional  shape  and  the 
flow  conditions.  It  is  desirable  to  model  the  body  vortex  wake 
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by  means  of  a  rational  method  capable  of  considering  a  variety  of 
body  shapes  over  a  wide  range  of  incidence  angles  and  Mach  num¬ 
bers  up  to  the  critical  speed.  It  is  important  that  the  separa¬ 
tion  vortex  wake  induced  effects  on  the  nonlinear  aerodynamic 
characteristics  of  the  body  be  handled  properly  with  a  method 
which  correctly  represents  the  physical  characteristics  of  the 
flow  field. 

The  phenomena  of  interest  are  the  sheets  of  vorticity  formed 
on  the  lee  side  of  the  body  at  high  angles  of  attack.  The  vorti¬ 
city  is  formed  by  boundary-layer  fluid  leaving  the  body  surface 
from  separation  points  on  both  sides  of  the  body  (Fig.  1).  At 
moderate  angles  of  incidence,  the  vorticity  rolls  up  into  a  sym¬ 
metrical  vortex  pair,  but  at  higher  angles,  the  vorticity  becomes 
asymmetric  as  shown  in  the  sketch  in  Figure  1(b).  A  method  to 
predict  these  flow  phenomena  in  the  vicinity  of  circular  and 
noncircular  bodies  in  subsonic  flow  is  described  in  References  1 
and  2.  The  extension  of  the  method  to  supersonic  flow  is  de¬ 
scribed  in  Reference  3. 

The  purpose  of  this  report  is  to  describe  an  engineering 
prediction  method  and  associated  computer  code  developed  to  cal¬ 
culate  the  nonlinear  aerodynamic  characteristics  and  flow  fields 
of  noncircular  bodies  at  high  angles  of  attack  at  speeds  up  to 
the  critical  speed.  The  objectives  of  the  method  are  to  use  a 
three-dimensional,  attached  flow,  potential  method  to  represent 
the  body  and  a  two-dimensional,  incompressible,  separated  flow 
model  to  calculate  the  lee  side  vortex  shedding  from  the  body 
alone  at  angle  of  attack  and  angle  of  roll.  The  predicted  pres¬ 
sure  distribution  on  the  body  under  the  influence  of  the  free 
stream  and  the  separation  vortex  wake  is  used  to  calculate  the 
aerodynamic  loads  on  the  body.  Conformal  mapping  technigues  are 
used  to  transform  noncircular  cross  sections  to  a  circle  for 
calculation  purposes. 
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The  following  sections  of  this  report  include  a  discussion 
of  the  approach  to  the  problem  and  a  description  of  the  analysis 
and  flow  models  required  to  carry  out  the  calculation.  The  pre¬ 
diction  method  is  evaluated  through  comparison  of  measured  and 
predicted  results  on  a  variety  of  body  shapes,  including  circular 
and  elliptical  cross  sections.  A  user's  manual  for  the  computer 
code  is  also  included.  The  manual  consists  of  descriptions  of 
input  and  output  and  of  sample  cases. 
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surface  velocity  in  crossflow  plane 
axial  perturbation  velocity 
local  velocity 

radial  perturbation  velocity 
vortex  induced  velocity 
velocity  components  in  real  plane 
free-stream  velocity 
complex  potential.  Eg.  (26) 
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exterior  angles  of  body  segment  for  numerical  map¬ 
ping,  Eq.  (6) 

vortex  strength 

axial  length  increment 

complex  coord i ~  te  in  an  intermediate  plane. 

Fig.  2(b) 

vertical  coordinate  in  an  intermediate  plane, 

Fig.  2(b) 

polar  angle  in  v-  plane.  Fig.  2(a) 

complex  coordinate  in  circle  plane.  Fig.  2(a);  also 
kinematic  viscosity 

run  length,  Eqs.  (46)  and  (47);  or  lateral  coordi¬ 
nate  in  an  intermediate  plane.  Fig.  2(b) 

free-stream  density 

complex  coordinate  in  real  plane.  Fig.  2(a) 

lateral  and  vertical  coordinates  in  circle  plane. 
Fig.  2(a) 

roll  angle 

velocity  potential  in  real  plane 


stream  function  in  real  plane 
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center  of  pressure 
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GENERAL  APPROACH 


Bodies  at  high  angles  of  attack  exhibit  distributed  vorti- 
city  fields  on  their  lee  side  due  to  boundary-layer  fluid  leaving 
the  body  surface  at  separation  lines.  One  approach  for  modeling 
these  distributed  vorticity  fields  has  involved  the  use  of  clouds 
of  discrete  potential  vortices.  Underlying  the  basic  approach  is 
the  analogy  between  two-dimensional  unsteady  flow  past  a  body  and 
the  steady  three-dimensional  flow  past  an  inclined  body.  In 
fact,  the  three-dimensional  steady  flow  problem  is  reduced  to  the 
two-dimensional  unsteady  separated  flow  problem  for  solution. 
Linear  theory  for  the  attached  flow  model  and  slender  body  theory 
to  represent  the  interactions  of  the  vortices  are  combined  to 
produce  a  nonlinear  prediction  method.  The  details  of  the  appli¬ 
cation  of  this  approach  to  the  prediction  of  subsonic  flow  about 
circular  and  noncircular  bodies  are  presented  in  References  1  and 
2.  Other  investigators  have  used  this  approach  to  successfully 
model  the  subsonic  flow  phenomena  in  the  vicinity  of  circular 
cross  section  bodies  (Refs.  4  and  5).  The  purpose  of  this  report 
is  to  document  in  detail  the  subsonic  analysis  of  References  1 
and  2  and  to  describe  the  code  VTXCLD  assembled  to  accomplish  the 
subsonic  vortex  shedding  calculations. 

The  calculation  procedure  is  carried  out  in  the  following 
manner.  The  axisymmetric  body  is  represented  by  discrete  point 
sources  and  doublets,  and  the  strength  of  the  individual  singu¬ 
larities  is  determined  to  satisfy  a  flow  tangency  condition  on 
the  body  in  a  nonseparated  uniform  flow  at  angle  of  incidence  and 
roll.  Compressibility  effects  on  the  body  are  included  by  a 
Gothert  transformation  which  keeps  the  cross  section  shape  un¬ 
changed  but  stretches  the  axial  body  coordinate.  Starting  at  a 
crossflow  plane  near  the  body  nose,  the  pressure  distribution  on 
the  body  is  computed  using  the  full  compressible  Bernoulli  egua- 
tion.  The  boundary  layer  in  the  crossflow  plane  is  examined  for 
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separation  using  modified  versions  of  Stratford's  separation 
criteria.  The  Stratford  separation  criteria  are  based  on  two- 
dimensional  incompressible  flow.  At  the  predicted  separation 
points,  incompressible  vortices  with  their  strength  determined  by 
the  vorticity  transport  in  the  boundary  layer  are  shed  into  the 
flow  field.  The  trajectories  of  these  free  vortices  between  this 
crossflow  plane  and  the  next  plane  downstream  are  calculated  by 
integration  of  the  equations  of  motion  of  each  vortex,  including 
the  influence  of  the  free  stream,  the  body,  and  other  vortices. 
The  pressure  and  trajectory  calculations  are  carried  out  by  map¬ 
ping  the  noncircular  cross  section  shape  to  a  circle  using  either 
analytical  or  numerical  conformal  transformations.  The  vortex 
induced  velocity  contribution  to  the  body  tangency  boundary  con¬ 
dition  includes  image  vortices  in  the  circle  plane.  At  the  next 
downstream  crossflow  plane,  new  vortices  are  shed,  adding  to  the 
vortex  cloud  representing  the  wake  on  the  lee  side  of  the  body. 
This  procedure  is  carried  out  in  a  stepwise  fashion  over  the 
entire  length  of  the  body.  The  details  of  the  individual  methods 
combined  into  the  prediction  method  are  described  in  the  fol¬ 
lowing  section. 


METHODS  OF  ANALYSIS 

The  development  of  an  engineering  method  to  predict  the 
pressure  distributions  on  arbitrary  missile  bodies  in  subsonic 
flow  at  high  angles  of  incidence  requires  the  joining  of  several 
individual  prediction  methods.  In  this  section,  the  individual 
methods  are  described  briefly,  and  the  section  concludes  with  a 
description  of  the  complete  calculation  procedure. 

Conformal  Mapping 

The  crossflow  plane  approach  applied  to  arbitrary  missile 
bodies  results  in  a  noncircular  cross  section  shape  in  the  pres- 
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ence  of  a  uniform  crossflow  velocity  and  free  vortices  in  each 
plane  normal  to  the  body  axis.  The  procedure  used  to  handle  the 
noncircular  shapes  is  to  determine  a  conformal  transformation  for 
mapping  every  point  on  or  outside  the  arbitrary  body  to  a  corre¬ 
sponding  point  on  or  outside  a  circular  body.  The  two- 
dimensional  potential  flow  solution  around  a  circular  shape  in 
the  presence  of  a  uniform  flow  and  external  vortices  is  well 
known  and  has  been  documented  numerous  places  in  the  literature 
(Refs.  6  and  7).  Thus,  the  procedure  is  to  obtain  the  potential 
solution  for  the  circular  body  and  transform  it  to  the  noncir¬ 
cular  body.  Conformal  transformations  used  are  of  two  distinct 
types,  analytical  and  numerical. 

Analytic  transformation.-  For  very  simple  shapes  like  an 
ellipse  [Fig.  2(a)],  the  transformation  to  the  circle  can  be 
carried  out  analytically  as  described  in  Reference  7.  For 
example. 


where 


a  =  v  + 


-  b' 


4  v 


(1) 


a  =  y  +  iz 


in  the  real  plane,  and 


v  =  t  +  i  X 


(2) 


(3) 


in  the  circle  plane. 


The  derivative  of  the  transformation  is 


(4) 
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which  is  required  for  the  velocity  transformation  discussed 
later. 


Numerical  transformation.-  For  complex  noncircular  shapes, 
the  transformation  cannot  be  carried  out  analytically  and  a  nu¬ 
merical  transformation  is  required.  The  numerical  transformation 
chosen  is  described  in  detail  in  Reference  8,  and  a  brief  summary 
of  the  conformal  mapping  procedure  follows. 

The  sequence  of  events  in  the  numerical  mapping  is  shown  in 
Figure  2(b).  The  arbitrary  cross  section  shape  of  the  body  in 
the  o-plane  is  required  to  have  a  vertical  plane  of  symmetry. 
The  transformation  of  interest  will  map  the  region  on  and  outside 
the  body  in  the  o-plane  to  the  region  on  and  outside  a  circle  in 
the  v-plane.  The  first  step  is  a  rotation  to  the  c-plane  so 
that  the  cross  section  is  symmetric  about  the  real  axis.  Thus, 

C  =  io 

=  5  +  in  (5) 

A  mapping  that  transforms  the  outside  of  the  body  in 
the  c-plane  to  the  outside  of  the  unit  circle  is 

Tr 

.  m  - 

C  =  /  — 2  n(v_vr)ir  dv  (6) 

v  r=l 


where  yr  are  the  exterior  angles  of  the  m  segments  of  the  body 
cross  section.  In  Equation  (6),  n  denotes  a  product  series.  For 
a  closed  body 

m 

l  Yr  «  2n  (7) 

r=l 
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The  final  transformation  has  the  form 


o 


-i 


(8) 


where  the  A  coefficients  are  obtained  through  an  iterative 

X> 

scheme  described  in  Reference  8  and  rQ  is  the  radius  of  the 
equivalent  circle  in  the  v-plane.  The  derivative  of  the  trans¬ 
formation  is 


do 
d  v 


-i 


1 


L 


I 

4*1 


r£+1  l A 
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l+l 


(9) 


which  is  required  for  the  velocity  calculations  described  in  a 
later  section. 

The  above  numerical  mapping  procedure  has  been  applied  to  a 
wide  range  of  general  cross  section  shapes  with  good  success. 

Body  Model 

A  three-dimensional  representation  of  the  missile  volume  is 
needed  for  purposes  of  predicting  the  absolute  pressure  coeffi¬ 
cient  on  the  surface  (Ref.  2).  Since  the  model  must  be  a  compu¬ 
tationally  efficient  means  of  representing  both  circular  and 
noncircular  bodies  in  compressible  flow  up  to  the  critical  speed, 
a  method  using  discrete  singularities  on  the  body  axis  is  de¬ 
scribed  in  this  section.  As  described  in  Reference  3,  a  panel 
method  can  be  used  to  represent  the  missile  surface;  however, 
increased  computational  requirements  make  the  present  approach 
more  desirable  for  an  engineering  prediction  method. 

Circular  bodies.-  The  volume  of  an  axisymmetric  body  is  well 
represented  by  a  series  of  point  sources  and  sinks  distributed  on 
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the  axis.  A  number  of  models  in  varying  degrees  of  complexity 
are  available  for  this  task;  for  example,  see  References  9-12. 
For  use  in  VTXCLD,  a  discrete  source  model  from  the  latter  refer¬ 
ence  was  selected  for  its  accuracy,  economy,  and  reliabilty. 


Given  a  series  of  K  point  sources  (sinks)  distributed  on  the 
missile  axis,  the  induced  axial  and  radial  velocities  at  a  point 
( x ,  r )  are 


where 


(10) 


(11) 


(12) 


is  the  dimensionless  source  strength.  0£  and  xk  are  the  source 
strength  and  axial  locations,  respectively,  of  the  k-th  source. 
The  body  surface  slope  at  the  j-th  body  point  obtained  from  Equa¬ 
tions  (10)  and  (11)  is 


dr 

dx 


j 


h. 


V  cos  a 

or  c 


U. 


1  + 


V  cos a. 


for  j=l 


( N-3 ) 


(13) 
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Equation  (13)  produces  a  set  of  K-3  linear  equations  in  the  un¬ 
known  source  strengths.  For  a  closed  body,  imposing  the  condi¬ 
tion  of  the  sum  of  the  source  strengths  to  be  zero  gives  the 
relation 

K 

l  0/  =  0  (14) 

k  =  l  K 

The  remaining  two  conditions  are  the  imposition  of  stagnation 
points  at  the  nose  and  tail  of  the  body.  These  conditions  from 
Equation  (10)  are 


Given  the  source  positions  on  the  missile  axis.  Equations  (13) 
through  (16)  comprise  the  total  set  of  equations  to  solve  for  the 
source  strengths.  The  predicted  hull  shape  from  the  stream  func¬ 
tion  is 


which  must  be  solved  by  iteration. 

The  above  method  has  proved  successful  in  modeling  a  variety 
of  axisymmetric  missile  bodies;  however,  some  care  is  required 


because  of  the  ill-conditioned  matrix.  Best  results  are  obtained 
if  the  discrete  sources  are  spaced  at  intervals  of  60%  of  the 
local  radius  and  if  the  surface  slope  description  of  the  body  is 
smooth  without  discontinuities.  Boundary  conditions  are  satis¬ 
fied  at  points  midway  between  the  source  locations.  The  predic¬ 
tion  of  an  approriate  source  distribution  is  an  automated  part  of 
VTXCLD,  and  the  user  must  only  specify  missile  geometry. 

Noncircular  bodies.-  An  appropriate  body  model  for  missiles 
with  noncircular  cross  sections  is  a  surface  panel  method  similar 
to  that  described  in  Reference  3;  however,  the  use  of  a  panel 
model  adds  significantly  to  the  cost  of  each  computation.  For 
this  reason  an  alternate  approach  was  selected  for  use  with  non¬ 
circular  bodies. 

The  noncircular  body  is  replaced  with  an  equivalent  axisym- 
metric  body  having  the  same  cross  sectional  area  distribution  as 
the  actual  body.  There  are  approximations  involved  with  this 
model  since  the  induced  u-velocity  due  to  noncircular  effects  are 
neglected;  however,  for  bodies  that  do  not  deviate  greatly  from 
circular,  the  approximations  are  acceptable.  Some  guidelines  are 
included  with  the  user's  instructions  to  illustrate  possible 
problem  areas  and  the  magnitude  of  the  approximation  associated 
with  the  model. 

Compressibility  effects.-  The  selection  of  a  compressibility 
correction  scheme  must  take  into  consideration  the  configurations 
of  interest  and  the  calculation  procedure.  For  example,  vortex 
shedding  and  tracking  is  very  dependent  on  cross  sectional  shape, 
and  numerical  transformations  add  significantly  to  computation 
time;  therefore,  it  is  advantageous  that  the  compressibility 
transformation  used  have  no  effect  on  the  cross  sectional 
shape.  With  this  guideline  and  based  on  similar  experiences  and 
requirements  (Ref.  12),  a  transformation  which  modifies  only  the 
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axial  coordinate  is  needed.  The  Gothert  Rule,  described  in  Ref¬ 
erence  13  and  others,  is  the  choice  for  the  compressibility 
transformation.  A  brief  description  of  the  method  included  in 
VTXCLD  follows. 

The  transformation  from  the  compressible  (x,y,z)  coordinate 

system  to  the  incompressible  (x',y',z')  system  is 

x'  =  x/3 

y'  =  y 

z'  =  z  (18) 


where 


3  -  /i-m2  (19) 

or 

As  a  result  of  Equation  (18),  the  body  slope  and  angle  of  attack 
become 


dr'  -  8 

dx'  3 

dr 

(20) 

dx 

a'  =  3a 

(21) 

and  the  velocity  fields  are  related  by 

u  =  u’/B2 
v  =  v'/B 

w  =  w' /B  ( 22 ) 
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The  vortex  shedding  scheme  in  VTXCLD  requires  that  the  compressi¬ 
bility  correction  be  applied  in  a  manner  sliqhtly  different  from 
the  usual.  At  the  initiation  of  the  calculation,  the  actual 

missile  body  is  transformed  to  the  incompressible  shape  by  the 

above  stretching  procedure,  and  the  source  distribution  is  ob¬ 

tained  as  described  above.  The  modified  incompressible  flow 
conditions  are  used  to  calculate  the  velocity  field  associated 
with  the  body.  This  velocity  field  is  transformed  back  to  the 

compressible  domain  so  that  it  is  available  for  use  in  vortex 
tracking  and  compressible  pressure  coefficient  calculations.  The 
compressible  pressure  coefficient  is  used  to  locate  separation 
and  define  the  shed  vorticity.  Since  the  compressibility  trans¬ 
formation  has  no  direct  effect  on  the  cross  sectional  shape  or 
the  discrete  wake  vortices,  the  separation  and  vortex  tracking 
calculations  take  place  in  the  compressible  domain.  No  compres¬ 
sibility  effects  are  considered  for  the  vortex  induced  velo¬ 
cities  . 


Vortex  Shedding  Model 

The  body  vortex  shedding  model  described  in  this  section  is 
nearly  the  same  as  the  subsonic  model  presented  in  Reference  2; 
therefore,  only  the  differences  between  the  two  models  will  be 
discussed  in  detail. 


Equations  of  motion.-  The  equations  o*  motion  of  a  shed  nose 
vortex  in  the  presence  of  other  free  vortices  in  the  vicinity  of 
a  body  in  a  uniform  stream  follow.  In  the  circle  (v)  plane,  the 
position  of  a  vortex,  r  ,  is 


V 

m 


r 


T  + 

m 


and  an  imaqe  of  r  is  located  at 

m 


(23) 
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(24) 


-r  v 


to  satisfy  the  flow  tangency  condition  on  the  body  surface. 
In  the  real  plane,  the  position  of  the  vortex  Tm  is 


a  =  y  +  iz 

ntj,  m 

The  complex  potential  in  the  real  plane  is 


(25) 


W{  a)  =  $  —  i 


(26) 


and  the  corresponding  velocity  at  Fm  is 

dVo)  d  d“l 

vm  ‘  lwm  *  '  33  ~  =  35  Wm  °  33 1  o=. 


(27) 


The  complex  potential  of  is  not  included  in  Eguation  (27)  to 
avoid  the  singularity  at  that  point.  The  derivative  of  the 
transformation  is  obtained  from  Eguations  (4)  or  (9). 

The  total  velocity  at  r  in  the  crossflow  plane  is  written 

m 


v  -  xw 
it  m 


- - -  =  G  +  G  +  G  +  G  +  G_  +  G 

V  a  p  n  m  T  r 

or 


(28) 


where  each  term  in  Eguation  (28)  represents  a  specific  velocity 
component  in  the  o-plane.  The  first  term  represents  the  uniform 
flow  due  to  angle  of  attack 


G  =  -i  sina  1  + 
a 


rlol  dv| 

LVJ  d0'a=ar 


(29) 
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The  second  term  represents  the  uniform  flow  due  to  angle  of  yaw. 


G0  *  -st 


ins  1  -  [^1 

|v  I  da ' a=a 
L  m  J  r 


(30) 


In  compressible  flow,  the  velocity  components  from  Equations  (29) 
and  (30)  used  in  the  surface  pressure  calculations  include  ef¬ 
fects  of  the  compressibility  transformation  described  in  a  pre¬ 
vious  section. 


The  next  term  represents  the  influence  of  all  vortices  and 
their  images,  with  the  exception  of  r  . 

N  rn  x 

n=l  2*roV«  (vm/ro)-(ro/vn) 


(31) 


(v  /r  )  —  ( v  /r  )  1  da 1 

mo  no  n*m  a-a. 


The  next  term  is  due  to  the  image  of  r  . 

m 


G  =  i  3-V  - - - —  ^1 

m  *ro  ~  ( v  /c  ) - ( n  / v  )  °  0=0 

m  o  o  m  m 


(32) 


The  fifth  term  in  Equation  (28)  represents  the  potential 


of  r  in  the  a-plane  and  is  written  as 
m 


r  =  -i  -I*-  m  m 

T  2w V  L2J  dv  LdcJ 


(33) 


The  last  term  in  Equation  (28)  represents  the  velocity  components 
induced  by  the  body  source  singularities  representing  the  volume 
effects . 
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(34) 


where  vr  and  wr  are  components  of  vR  from  Equation  (11). 


Since  the  body  source  singularities  are  three-dimensional, 
they  contribute  an  induced  axial  velocity,  ur,  given  by  Equation 
(10).  Compressibility  effects  are  included  through  the  stretch¬ 
ing  of  the  axial  coordinate  and  subsequent  transformation  of  the 
induced  velocity  components. 


The  differential  equations  of  motion  for  r  are 

m 


v  -  1W 
m _ m 

V  cosa  +  u^. 
a-  c  r 


(35) 


where 


(36) 


Therefore,  the  two  equations  which  must  be  integrated  along  the 
body  length  to  determine  the  trajectory  of  are 


dx  V  cosa  +  u 
op  c  r 


(37) 


and 


dx  V  cosa  +  u 
op  c  r 


(38) 


There  are  a  pair  of  equations  like  (37)  and  (38)  for  each  vortex 
in  the  field.  As  new  vortices  are  shed,  the  total  number  of 
equations  to  solve  increases  by  two  for  each  added  vortex.  These 
differential  equations  are  solved  numerically  using  a  method 
which  automatically  adjusts  the  step  size  to  provide  the  speci¬ 
fied  accuracy. 
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Surface  pressure  distribution.-  The  surface  pressure  dis¬ 
tribution  on  the  body  is  required  to  calculate  the  forces  on  the 
body  and  the  separation  points.  The  surface  pressure  coefficient 
is  determined  from  the  Bernoulli  equation  in  the  form 


where 


and 


(39) 


(40) 


(41) 


where  U  is  the  total  velocity  (including  V  )  at  a  point  on  the 
body. 


The  last  term  in  Equation  (41)  represents  the  axial  velocity 
missing  from  the  two-dimensional  singularities  in  the  flow 
model.  In  this  case,  the  shed  vortices  and  the  doublet  part  of 
Equations  (29)  and  (30)  are  the  only  singularities  contributing 
to  this  term.  The  complex  potential  of  these  two-dimensional 
singularities  is 

W(v)  =  (j>  +  i  i|< 
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Separated  wake.-  The  separated  wake  on  the  lee  side  of  the 
body  is  made  up  of  a  large  number  of  discrete  vortices,  each 
vortex  originating  from  one  of  two  possible  locations  at  each 
time  step  in  the  calculation.  The  major  portion  of  the  lee  side 
vortex  wake  has  its  origin  at  the  primary  separation  points  on 
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each  side  of  the  body.  The  remainder  of  the  wake  originates  from 
the  secondary  separation  points  located  in  the  reverse  flow  re¬ 
gion  cn  the  lee  side  of  the  body.  Both  of  these  points  are  il¬ 
lustrated  in  the  sketch  of  a  typical  crossflow  plane  of  an  ellip¬ 
tic  cross  section  body  shown  in  Figure  3.  The  mechanics  of  the 
calculation  of  the  individual  vortices  follows. 

As  described  in  References  1  and  2,  the  pressure  distribu¬ 
tion  for  the  primary  flow  in  the  crossflow  plane  is  referenced  to 
the  conditions  at  the  minimum  pressure  point,  and  a  virtual  ori¬ 
gin  for  the  beginning  of  the  boundary  layer  is  calculated.  The 
adverse  pressure  distribution  downstream  of  the  minimum  pressure 
point  is  considered  with  either  Stratford's  laminar  (Ref.  14)  or 
turbulent  (Ref.  15)  separation  criterion  to  determine  whether  or 
not  separation  has  occurred.  These  criteria,  based  on  two- 
dimensional,  incompressible,  flat  plate  data,  are  adjusted  for 
three-dimensional  crossflow  effects  in  Reference  2.  The  laminar 
separation  criterion  states  that  the  laminar  boundary  layer  sepa¬ 
rates  when 


In  a  turbulent  boundary  layer,  separation  occurs  when 


(46) 


(Re^x  10 


-0.1 


0.35  sina 

c 


(47) 


The  sinac  in  Equations  (46)  and  (47)  is  the  three-dimensional 
modification. 

If  the  criteria  indicate  a  separation  point,  the  vorticity 
flux  across  the  boundary  layer  at  separation  is  shed  into  a 
single  point  vortex  whose  strength  is 
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(48) 


V  -  2  cos  a 

op  2V  c 


OP 

assuming  no  slip  at  the  wall.  6  is  the  vortex  reduction  factor 
described  below. 

The  initial  position  of  the  shed  vortex  is  determined  such 
that  the  surface  velocity  in  the  crossflow  plane  at  the  separa¬ 
tion  point  is  exactly  canceled  by  the  shed  vortex  and  its 
image.  When  this  criterion  results  in  a  vortex  initial  position 
that  is  too  near  to  the  body  surface,  certain  numerical  problems 
can  cause  difficulty  in  calculating  the  trajectory  of  this  vor¬ 
tex.  If  the  initial  position  of  the  vortex  is  nearer  than  five 
percent  of  the  body  radius  from  the  surface,  the  vortex  is  gener¬ 
ally  placed  five  percent  of  the  equivalent  body  radius  from  the 
surface.  This  position  represents  the  approximate  thickness  of 
the  boundary  layer. 

The  calculation  of  secondary  separation  is  carried  out  in 
the  same  manner  as  described  above  for  primary  separation.  It  is 
necessary  that  a  reverse  flow  region  exist  on  the  lee  side  of  the 
body  and  that  a  second  minimum  pressure  point  be  found  in  this 
region.  Surface  flow  visualization  of  secondary  separation  re¬ 
gions  on  bodies  at  high  angles  of  attack  in  subsonic  flow  are 
shown  in  Reference  16.  For  purposes  of  this  analysis,  the  re¬ 
verse  flow  is  assumed  to  be  laminar  from  the  lee  side  stagnation 
point  to  the  secondary  separation  point,  and  Stratford's  laminar 
criterion  is  used  to  locate  secondary  separation.  Laminar  sepa¬ 
ration  in  the  reverse  flow  region  is  expected  because  of  the  low 
velocities  on  this  portion  of  the  body.  The  vortex  released  into 
the  flow  at  the  secondary  separation  point  has  the  opposite  sign 
of  the  primary  vortex  and  is  generally  weaker  in  strength. 

Transition .-  In  general,  the  primary  separation  on  a  missile 
is  either  laminar  or  turbulent,  depending  on  the  Reynolds  number 
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of  the  flow;  however,  as  noted  in  Reference  17,  there  is  an  angle 
of  attack  and  Reynolds  number  range  over  which  the  character  of 
the  separation  is  transitional.  Flow  visualization  of  the  tran¬ 
sition  phenomenon  on  missile  shapes  is  available  in  References  16 
and  18.  Body  vortex  shedding  technology  is  not  developed  to  the 
level  that  prediction  of  the  beginning  of  transition  nor  the 
length  of  the  transition  is  practical.  Correlation  of  the  lim¬ 
ited  amount  of  data  available  has  permitted  the  authors  to  define 
approximately  the  limits  of  transition.  The  results  of  the  cor¬ 
relation  indicate  that  the  separation  switches  from  laminar  to 
turbulent  when  the  Reynolds  number  based  on  the  length  from  the 
nose  times  the  sine  of  the  angle  of  incidence  is  approximately 
7x10^.  The  length  of  the  transition  region  is  about  one  body 
diameter.  More  analysis  is  reguired  to  verify  the  correlation, 
and  the  authors  recommend  that  transition  effects  be  included 
with  great  care. 

Reduction  factor.-  The  empirical  vortex  reduction  factor, 

6,  in  Equation  (48)  is  widely  used  in  discrete  vortex  methods  to 
provide  better  agreement  between  measured  and  predicted  vorticity 
strength  in  the  wake.  Various  investigators  (e.g..  Refs.  2,  4, 
and  19)  recommend  a  value  in  the  range  0.6  <  6  <  1.0  for  subsonic 
flow  with  0.6  being  the  more  common  value.  Based  on  the  experi¬ 
ence  of  the  authors,  the  proper  value  of  6  can  vary  with  the 
geometry  and  flow  conditions.  For  example,  6  =  0.6  is  an  appro¬ 
priate  value  for  an  ogive-cylinder  missile  in  subsonic  flow, 
but  6  =  1.0  is  a  better  value  for  a  boattail  body  in  subsonic 
flow.  Unless  otherwise  noted,  6  =  0.6  is  used  for  all  predicted 
results  presented  in  this  report. 

Forces  and  moments.-  The  forces  and  moments  on  the  body  are 
computed  by  integration  of  the  pressure  distribution  around  the 
body.  At  a  specified  station  on  the  body,  the  normal-force  coef¬ 
ficient  on  a  Ax  length  of  the  body  is 
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where  r'  is  the  distance  from  the  axis  of  the  body  to  the  body 
surface  and  g'  is  the  local  slope  of  the  body  in  the  crossflow 
plane.  This  is  illustrated  in  the  sketch  in  Figure  4.  For  cir¬ 
cular  bodies,  r'  =  rQ  and  6'  =3.  The  total  normal  force  coeffi¬ 
cient  on  the  body  is 
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and  the  pitching-moment  coefficient  is 
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The  center  of  pressure  of  the  normal  force  is 
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Similarly,  the  side-force  coefficient  on  a  Ax  length  of  the 
body  is 


[AS]  2ir 

'  4d  -  -  5  /  V  sin8' 


dg 


The  total  side-force  coefficient  on  the  body  is 

f  r 
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and  the  ya wing-moment  coefficient  is 
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The  center  of  pressure  of  the  side  force  is 
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(56) 


A  noncircular  body  at  an  arbitrary  roll  angle  may  experience 
a  rolling  moment  caused  by  the  nonsymmetry  of  the  loading  around 
the  body.  The  total  rolling  moment  on  the  body  is  calculated  by 
summing  the  moments  of  the  individual  components  of  normal  force 
and  side  force  around  each  cross  section  and  integrating  over  the 
body  length.  Using  Equations  (4)  and  (5),  the  rolling  moment 
coefficient  is  calculated  as 
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Though  the  primary  purpose  of  the  vortex  shedding  method  is 
not  to  predict  drag  or  axial  force  on  the  missile,  a  procedure  to 
estimate  both  pressure  drag  (excluding  base  drag)  and  skin  fric¬ 
tion  is  included.  The  pressure  contribution  to  the  axial-force 
coefficient  is 
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The  local  skin-friction  coefficient,  based  on  the  assumption  of  a 
l/7th  power  law  velocity  profile  in  the  boundary  layer,  is 
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which  produces  a  drag  coefficient  due  to  friction 
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At  high  angles  of  incidence,  streamlines  around  the  body  are 

inclined  at  approximately 

o  =  tan^  (2  tan  a  )  (61) 

s  c 

therefore,  the  axial  component  of  the  friction  drag  is 

C.  =  Cp.  cos  a  (62) 

A  j  Dg  s 

The  total  axial-force  coefficient  is  the  sum  of  Equations  (58) 
and  (62). 


Vortex  core.  -  The  diffusion  core  model  (Ref.  2)  for  the 
point-vortex  induced  velocities  removes  the  singularity  at  the 
vortex  origin  and  effectively  reduces  the  velocities  near  the 
vortex.  The  tangential  velocity  induced  by  a  single  point  vortex 


where  r  is  the  distance  from  the  vortex  to  the  field  point  and  x 
is  a  measure  of  the  age  of  the  vortex.  The  induced  velocity  from 
this  core  model  is  illustrated  in  the  following  sketch. 
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The  vortex  core  model  represented  by  Equation  (63)  has  re¬ 
ceived  considerable  attention  in  the  context  of  body  vortex  in¬ 
duced  effects,  and  it  has  a  number  of  shortcomings.  Since  the 
exponential  term  is  a  function  of  r,  the  flow  medium  (v),  and  the 
age  of  the  vortex  (x),  the  core  is  constantly  changing  size  as 
the  vortex  moves  through  the  field.  Under  certain  conditions, 
the  core  radius,  denoted  as  rc  in  the  sketch,  can  become  very 
small  and  the  induced  velocity  becomes  unrealistically  large.  In 
an  attempt  to  further  modify  the  core  model  to  keep  the  induced 
velocities  within  physically  realistic  limits,  the  following 
modification  was  made  in  Reference  3.  The  location  of  the  maxi¬ 
mum  induced  velocity,  rc  in  the  sketch,  is  fixed  at  a  specific 
radius  to  be  selected  by  the  user.  Given  a  core  radius,  the 
vortex  induced  velocity  is 


The  core  model  described  by  Equation  (64)  is  included  in  the 
version  of  the  code  described  in  this  report,  and  in  a  later 
section,  guidelines  for  the  selection  of  an  appropriate  core 
radius  are  presented.  Results  indicate  that  this  simple  core 
model  provides  adequate  smoothing  for  the  necessary  velocity 
calculations. 

Other  investigations  of  discrete  vortex  models  must  use  a 
core  model  of  some  type.  Although  there  are  many  different  core 
models,  all  serve  the  same  purpose  as  those  described  above,  and 
nearly  all  are  directed  at  eliminating  the  singularity  and  reduc¬ 
ing  the  maximum  induced  velocity.  The  code  VTXCLD  described  in 
this  report  is  easily  modified  to  incorporate  another  core  model 
if  the  user  desires  to  do  so. 

Asymmetry . -  An  axisymmetric  missile  at  large  angles  of  at¬ 
tack  in  subsonic  flow  will  exhibit  asymmetries  in  the  vortex 
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field  and  associated  induced  side  force.  The  vortex  cloud  ap¬ 
proach  described  herein  will  predict  no  asymmetry  at  any  angle  of 
attack  on  a  missile  with  a  vertical  plane  of  symmetry  unless  the 
solution  is  perturbed.  Every  investigator  (e.g.,  Refs.  2,  5,  and 
19)  has  a  favorite  means  of  introducing  a  perturbation  to  the 
solution,  and  no  good  argument  can  be  made  for  one  over  an¬ 
other.  Unfortunately,  the  size  and  type  of  perturbation  has  an 
effect  on  the  final  solution;  therefore,  some  empircal  informa¬ 
tion  is  reguired  to  calibrate  the  perturbation  for  a  particular 
geometry  and  flow  condition. 

In  Reference  2  and  in  VTXCLD,  the  perturbation  is  introduced 
by  a  small  modification  of  the  predicted  symmetric  separation 
points  on  either  side  of  the  missile  at  the  first  few  integration 
steps  near  the  nose.  The  solution  is  then  allowed  to  develop  on 
its  own.  Some  results  for  various  magnitude  perturbations  are 
presented  in  Reference  2  and  this  report. 

PROGRAM  VTXCLD 


Program  VTXCLD  is  a  rational  flow  model  which  predicts  the 
nonlinear  forces  and  moments  acting  on  slender  circular  and  non¬ 
circular  bodies  at  high  angles  of  attack  in  steady  subsonic  flow 
conditions.  The  nonlinear  effects  are  obtained  by  modeling  the 
lee  side  separation  vorticity  with  discrete  vortices  in  the 
crossflow  planes.  The  code  has  application  as  an  engineering 
prediction  method  for  use  in  preliminary  design  and  analysis 
studies  of  missiles  at  high  angles  of  incidence,  and  as  such,  it 
is  important  that  the  code  be  simple  and  economical  to  use  and 
reguire  a  minimum  of  empirical  information.  However,  it  is 
egually  important  that  the  code  adequately  represent  the  impor¬ 
tant  nonlinear  flow  phenomena  and  simulate  the  major  physical 
features  of  the  flow  field. 
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The  following  sections  are  a  user's  manual  for  program 
VTXCLD.  These  sections  include  descriptions  of  the  program,  the 
calculation  procedure,  input  preparation,  sample  input  and  output 
cases  which  illustrate  various  program  features,  and  program 
limitations. 


General  Description 

This  section  provides  a  general  description  of  program 
VTXCLD  and  its  various  subroutines.  The  code  is  written  in  modu¬ 
lar  form  so  that  as  improvements  in  flow  models  become  available 
or  as  other  modifications  are  required  they  are  easy  to  incorpo¬ 
rate.  The  flow  of  the  calculation  is  described  in  this  report 
and  detailed  comments  are  provided  throughout  the  code  to  assist 
the  user  in  understanding  the  order  of  calculation. 

The  computer  code  consists  of  the  program  VTXCLD  and  39 
subroutines.  The  overall  flow  map  of  the  program  is  shown  in 
Figure  5  where  the  general  relationship  between  the  subroutines 
and  external  references  can  be  seen.  Communication  between  the 
program  modules  is  handled  primarily  by  named  common  blocks.  A 
cross  reference  table  for  the  calling  relationship  between  the 
program  subroutines  is  shown  in  Figure  6.  A  similar  table  for 
the  named  common  blocks  is  shown  in  Figure  7. 

The  program  is  written  in  standard  FORTRAN.  Execution  on  a 
VAX  11/750  can  vary  from  10  seconds  to  more  than  45  minutes  de¬ 
pending  on  various  factors  such  as  the  geometry,  the  integration 
interval  (DX),  the  number  of  shed  vortices,  the  use  of  flow  sym¬ 
metry,  and  the  flow  conditions. 
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Subroutine  Description 


This  section  briefly  describes  the  main  program  VTXCLD  and 
its  various  subroutines. 

VTXCLD  This  is  the  main  program  for  calculating  the  forces 
and  moments  on  a  slender  body  in  a  steady  flow 
condition  including  the  nonlinear  effects  of  lee 
side  separation  vorticity. 

AMAP  Transforms  *  specified  point  in  the  circle  plane  to 

its  corresponding  point  in  the  body  plane. 

ASUM  Calculates  velocity  term  due  to  transformation  for 

noncircular  bodies. 

BMAP  Sets  up  a  table  of  points  on  the  actual  body  cor¬ 

responding  to  specified  points  on  the  circle. 

BODY  Organizes  the  flow  through  VTXCLD.  It  organizes 

calls  to  predict  vortex  shedding,  to  calculate 
surface  velocities  and  pressures,  and  to  calculate 
the  overall  forces  and  moments  on  the  body. 

CMAP  Transforms  a  specified  point  in  the  body  plane  to 

its  corresponding  point  in  the  circle  plane. 

COMBIN  Combines  vortices  which  are  separated  by  a  distance 
less  than  or  equal  to  "RGAM".  (The  resulting  vor¬ 
tex  strength  is  the  sum  of  the  individual  strengths 
and  its  location  is  the  centroid  of  the  combined 
vortices ) . 
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CON FOR 


CPPOT 


DFEQKM 


DPHIDT 


DZDNU 


F 


FDPSDR 


FORCE P 


FPSI 


FPVEL 


Calculates  the  coefficients  of  the  numerical  con¬ 
formal  mapping  function.  Transforms  the  region 
outside  of  a  polygonal  shape  with  a  vertical  plane 
of  symmetry  to  the  region  outside  of  a  circle. 

Calculates  a  potential  pressure  distribution  at  a 
given  x  station  (the  effects  of  separation  vortices 
are  not  present  in  the  Cp  calculation).  This  is 
used  in  determining  a  turbulent  separation  point 
for  the  transition  region. 

Obtains  vortex  paths  between  stations  x  and  x  to  Dx 
using  a  Kutta-Merson  integration  scheme. 

Calculates  the  two-dimensional  unsteady  pressure 
term  for  use  in  the  unsteady  Bernoulli  equation. 

Calculates  the  differential  of  the  transform  of  a 
polygonal  shape  to  a  circle. 

Calculates  the  derivatives  of  the  vortex  equations 
of  motion  which  are  used  in  conjunction  with  DFEQKM 
to  obtain  the  vortex  paths. 

Calculates  the  axial  velocity  due  to  a  three- 
dimensional  source  distribution. 

Calculates  the  forces  and  moments  on  the  body  by  an 
integration  of  the  surface  pressure  distribution. 

Calculates  the  body  shape  represented  by  a  three- 
dimensional  source  distribution. 

Calculates  the  velocity  at  specified  field  points. 
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INPT 

INVERS 

PLOT 

RDRDX 

RDRDXI 

SEPRAT 

SHAPE 

SKIN 

SORDIS 

SORVEL 


Reads  in  and  prints  out  the  body  geometry,  the  flow 
conditions,  various  program  options,  and  restart 
quantities . 

Solves  simultaneous  equations.  Ax  =  B. 

Creates  a  line  printer  plot  of  the  body  and  its 
vortex  wake  at  any  qiven  x  station.  It  is  composed 
of  the  followinq  subroutines:  PLOTV,  PL0TA2 , 
P LOTA 5 ,  PLOT A 6 ,  PL0TA7 ,  and  PL0TA8 . 

Calculates  r  and  dr/dx  at  a  qiven  x  station. 

Initializes  the  calculation  of  r  and  dr/dx  for 
subroutine  RDRDX. 

Predicts  the  position  in  the  crossflow  plane  where 
separation  occurs  based  on  the  Stratford  criteria 
for  laminar  and  turbulent  flows. 

Determines  the  body  geometric  parameters  from  a 
table  look-up  of  input  data. 

Estimates  the  axial  skin  friction  force  on  the 
body . 

Calculates  the  source  distribution  representing  an 
axisymmetric  body  of  revolution  given  the  body 
shape  and  slope  at  each  x  station. 

Calculates  perturbation  velocities  due  to  a  three- 
dimensional  source  distribution.  Used  by  subrou¬ 
tine  SORDIS  in  calculation  of  the  source  distribu¬ 
tion. 
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SRCVEL  Calculates  perturbation  velocities  due  to  a  three- 
dimensional  source  distribution. 

SUM  Analytical  mapping  of  an  elliptical  cross  section. 

TRANSN  Calculates  the  body  separation  points  in  the  tran¬ 

sition  region. 

VCENTR  Calculates  the  centroid  of  each  vortex  field. 

VLOCTY  Calculates  the  velocity  components  in  the  cross 
flow  plane  at  a  specified  field  point. 

VTABLE  Sets  up  the  arrays  of  free  vorticity  for  two 
cases:  (1)  The  vorticity  in  the  working  array  GAM 

is  distributed  into  the  individual  arrays  GAMP, 
GAMM,  GAMR,  GAMA,  or  (2)  the  vorticity  in  the  indi¬ 
vidual  arrays  GAMP,  GAMM,  GAMR,  GAMA  is  distributed 
into  the  working  array,  GAM,  in  the  proper  order. 

Z  Calculates  the  transformation  of  a  polygonal  shape 

to  a  circle  of  radius  RZ . 

Program  Limitations  and  Suggestions 

Program  VTXCLD  predicts  the  aerodynamic  forces  acting  on 
slender  bodies  alone  at  high  angles  of  attack,  with  some  limita¬ 
tions.  Certain  of  the  limitations  that  follow  are  basically 
rules-of-thumb  for  the  use  of  the  program  which  have  been  deter¬ 
mined  through  program  verification  and  experience  with  the  code. 

Mach  number . -  The  major  limitation  in  the  program  is  the 
Mach  number  or  compressibility  correction.  The  compressibility 
correction  added  to  VTXCLD  is  based  on  small  perturbation  theory 
which  restricts  the  body  to  be  slender,  the  angle  of  attack  to  be 
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small,  and  the  Mach  number  to  be  approximately  0.8  or  less.  The 
slender  body  restriction  is  not  a  new  restriction  as  this  assump¬ 
tion  is  inherent  in  the  crossflow  analogy  used  by  VTXCLD  to  pre¬ 
dict  lee  side  vorticity.  These  restrictions  result  from  the 
assumptions  made  for  small  perturbation  theory,  and  any  cases 
outside  these  restrictions  should  be  viewed  with  caution;  how¬ 
ever,  VTXCLD  has  been  applied  to  blunt  nose  shapes  such  as  spher¬ 
oids  or  ellipsoids  with  good  success.  Typical  missile  shapes 
with  slender  noses  pose  no  problems  to  the  method. 

Source  distribution.-  The  following  suggestions  should  be 
observed  when  using  the  source  distribution  calculation  in 
VTXCLD.  The  source  distribution  in  VTXCLD  places  three-dimen¬ 
sional  point  sources  on  the  axis  of  the  body  (or  equivalent  body 
in  a  noncircular  case).  For  subsonic  Mach  numbers  it  is  neces¬ 
sary  to  close  the  body  to  calculate  a  smooth  source  distribu¬ 
tion.  Generally,  the  body  can  be  closed  at  its  base  by  including 
a  smooth  boattail,  and  a  successful  technique  for  most  cases  is 
to  close  the  model  with  its  nose.  The  source  distribution  does 
not  handle  rapid  changes  or  discontinuities  in  body  slope  well. 

The  source  distribution  calculation  often  has  a  problem  with 
short  pointed  noses  as  well.  Short  blunted  noses  seem  to  cause 
no  difficulty  because  they  resemble  Rankine  bodies,  a  single 
point  source  in  a  uniform  stream.  Sharp  pointed  noses  often 
cause  an  oscillation  in  the  source  strength  on  the  afterbody 
resulting  in  an  oscillation  in  the  velocities  and  thus  an  oscil¬ 
lation  in  the  axial  pressure  distribution.  This  has  been  ob¬ 
served  with  a  two-caliber  ogive  nose  but  was  not  apparent  for  a 
2.5-caliber  ogive.  The  oscillations  were  minor  in  their  effect 
on  the  total  ;alculation  and  result,  but  the  user  should  be  aware 
of  their  possibility. 
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Incidence  angle.-  The  well-known  and  little-understood  phe¬ 
nomenon  of  wake  asymmetry  at  large  incidence  angles  (a  >  30°)  is 
available  as  an  option  in  VTXCLD.  As  described  in  a  previous 
section,  the  naturally  symmetric  nature  of  the  vortex  shedding 
analysis  must  be  perturbed  to  produce  an  asymmetric  solution; 
however,  the  asymmetric  capability  within  VTXCLD  has  been  demon¬ 
strated  and  will  be  described  in  the  RESULTS  section.  The  user 
must  be  aware  of  the  limitations  associated  with  this  portion  of 
the  code  in  order  to  apply  appropriate  judgement  in  the  evalua¬ 
tions  of  the  results.  Based  on  the  author's  experience  with 
asymmetry,  care  must  be  taken  in  the  choice  of  the  perturbation 
parameter  and  its  effect  on  the  results.  This  is  to  serve  as  a 
warning,  and  the  user  should  be  aware  that  this  portion  of  the 
code  is  in  need  of  additional  deve.lc.pment . 

Transition.-  The  transition  capability  included  in  the  code 
is  also  in  an  early  stage  of  development.  The  character  of  the 
vortex  wake,  whether  originating  from  laminar  or  turbulent  sepa¬ 
ration  points,  has  an  important  effect  on  the  predicted  nonlinear 
aerodynamics  characteristics.  As  previously  discussed,  addi¬ 
tional  development  is  needed  to  better  understand  the  origin  and 
length  of  the  transition  regions.  The  user  is  cautioned  to  be 
very  careful  in  the  interpretation  of  predicted  results  in  flow 
regimes  which  may  be  in  a  transition  region  as  defined  by  Lamont 
in  Reference  17. 

Secondary  separation.-  The  capability  to  predict  secondary 
separation  or  separation  in  the  reverse  flow  region  on  the  lee 
side  of  the  body  is  included  in  VTXCLD.  As  discussed  in  a  later 
section,  there  are  certain  flow  conditions  in  which  secondary 
separation  has  an  important  role;  for  example,  asymmetric  shed¬ 
ding  at  large  incidence  angles.  Numerical  problems  associated 
with  tracking  the  secondary  vorticity  have  been  observed,  and 
some  additional  effort  is  reguired  to  make  secondary  separation  a 
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useful  option.  Therefore,  it  is  recommended  that  this  option  be 
treated  as  a  research  feature  and  not  used  for  production  calcu¬ 
lations  . 


Error  Messages  and  Stops 

The  code  VTXCLD  has  numerous  internal  error  messages.  These 
messages  are  generally  explicit  and  are  described  below.  There 
are  several  execution  terminations  at  numbered  stops  within  the 
program.  These  are  described  in  the  table  below  with  some  sug¬ 
gestions  to  correct  the  problem. 

The  error  message  "CONVERGENCE  NOT  OBTAINED  IN  INTEGRATION" 
is  printed  out  when  the  trajectory  calculation  does  not  converge 
within  the  specified  tolerance,  E5 .  Also  printed  out  is  the 
axial  station,  XI,  and  the  integration  interval,  H.  This  error 
results  in  a  STOP  40  which  is  discussed  in  the  table  to  follow. 

In  the  calculation  of  the  source  distribution  it  is  possible 
to  get  the  message  "MATRIX  IS  SINGULAR."  This  occurs  when  one  of 
the  pivot  elements  in  the  influence  matrix  is  zero.  This  error 
results  in  a  STOP  300. 

The  other  major  error  messages  are  found  in  the  numerical 

.  ..  *  *  * 

conformal  mapping  routine  CONFOR.  The  message  ERROR  IN  THE 
SUM  OF  EXTERIOR  ANGLES***"  is  printed  if  the  routine  has  diffi¬ 
culty  calculating  appropriate  body  angles  within  a  specified 

tolerance.  Another  message  "***COEFFICIENT  AN(1)  =  _  GREATER 

THAN  +-  0.0001  AFTER  _  ITERATIONS  OUT  OF  _ "  is  printed  on 

the  rare  occasions  when  the  mapping  procedure  is  having  diffi¬ 
culty  converging  on  a  solution.  These  messages  are  informational 
and  do  not  result  in  a  program  STOP.  However,  the  user  is  ad¬ 
vised  to  check  the  transformed  shape  carefully  with  the  actual 
shape  before  using  the  results. 
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Termination 

Message 

STOP  1 

STOP  40 


STOP  41 

STOP  44 


Description 

This  is  the  normal  stopping  condition.  Pro¬ 
gram  VTXCLD 

The  vortex  trajectory  integration  has  failed 
to  converge  at  an  axial  station.  This  problem 
is  usually  caused  by:  (1)  two  or  more  vor¬ 
tices  in  close  proximity  rotating  about  one 
another  so  rapidly  that  the  integration  sub¬ 
routine  cannot  converge  on  a  solution,  or  (2) 
a  vortex  is  too  near  the  body  surface.  These 
situations  are  often  dependent  on  specific 
flow  conditions  and  the  body  shape.  The  fol¬ 
lowing  steps  should  be  tried  to  resolve  the 
problem.  increase  the  integration  tolerance 
E5,  or  increase  the  size  of  the  combination 
parameter  RGAM  so  that  it  is  greater  than  the 
distance  between  the  troublesome  vortices.  It 
often  helps  to  increase  the  integration  inter¬ 
val  DX  and  restart  the  calculation  procedure 
from  a  previous  station.  Subroutine  BODY 

The  vortex  trajectory  calculation  has  resulted 
in  a  vortex  inside  the  body.  This  is  a  rare 
but  serious  problem  which  indicates  a  major 
problem  with  the  solution.  Subroutine  BODY 

The  total  number  of  shed  vortices  exceeds 
200.  To  correct  this  combine  some  vortices 
and  restart  the  calculation,  or  increase  the 
integration  interval,  DX,  so  that  fewer  vor¬ 
tices  are  shed.  Subroutine  BODY 
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STOP 


STOP 


STOP 

STOP 


STOP 


STOP 


45  The  maximum  number  of  positive  or  negative 

primary  separation  vortices  exceeds  70.  Use 
the  same  remedy  as  for  STOP  44.  Subroutine 
BODY 

46  The  maximum  number  of  positive  or  negative 

secondary  (reverse  flow)  separation  vortices 
exceeds  30 .  Use  the  same  remedy  as  for  STOP 
44  but  only  combine  like  vortices.  If  one  of 
the  above  steps  does  not  work,  the  reverse 
flow  calculation  can  be  turned  off  by  setting 
NSEPR  =  0.  Subroutine  BODY 

201  There  is  an  error  in  the  vortex  tables  in 

202  subroutine  VTABLE.  Check  the  input  values  of 
NBLSEP ,  NSEPR,  NVP,  NVM ,  NVA ,  NVR,  and  RGAM  in 
Items  1  through  6. 

300  There  is  a  singular  influence  matrix  in  the 
solution  for  the  three-dimensional  source 
distribution.  Check  the  input  geometry  for 
errors.  Subroutine  INVERS 

301  The  maximum  number  of  iterations  is  exceeded 
in  the  source  calculation.  Check  the  input 
geometry.  This  problem  can  occur  if  the  body 
has  a  long  slender  nose  followed  by  a  cylin¬ 
drical  section  and  then  a  long  slender  tail. 
If  this  is  the  case,  one  solution  is  to  close 
the  body  off  with  a  blunter  and/or  shorter 
tail.  Subroutine  SORDIS 
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Input  Description 


This  section  describes  the  input  data  and  format  required  by 
program  VTXCLD.  The  input  formats  are  shown  in  Figure  8,  and  the 
definitions  of  the  individual  variables  are  provided  in  the  fol¬ 
lowing  paragraphs.  Note  that  some  length  parameters  in  the  input 
list  are  dimensional  variables;  therefore,  special  care  must  be 
taken  that  all  lengths  and  areas  are  input  in  a  consistent  set  of 
units. 

The  remainder  of  this  section  includes  a  description  of  each 
of  the  input  variables  and  indices  required  for  the  use  of  pro¬ 
gram  VTXCLD.  The  order  follows  that  shown  in  Figure  8  where  the 
input  format  for  each  item  and  the  location  of  each  variable  on 
each  card  is  presented.  Data  input  in  I-format  are  right  justi¬ 
fied  in  the  fields,  and  data  input  in  F-format  may  be  placed 
anywhere  in  the  field  and  must  include  a  decimal  point. 

Item  1  is  a  single  card  containing  15  integer  flags.  These 
flags  are  right  justified  in  a  five  column  field  and  are  used  to 
specify  several  options  available  in  VTXCLD. 


I  tern 


Variable 


Description 


1  (1515)  NCIR  Cross  section  flag,  determines  the 

type  of  bod"  cross  section. 

=  0,  circular  body 
=  1,  elliptical  body 
=  2,  arbitrary  body 


NCF  Mapping  flag,  specifies  whether  the 

required  mapping  is  input  or  calcu¬ 
lated.  (NCIR  =  2) 

=  0,  calculate  mapping 


-42- 


input  mapping 


ISYM 


NBLSEP 


NSEPR 


NDFUS 


NDPHI 


INP 


=  2, 

Vortex  shedding  symmetry  flag. 

=  0,  symmetric  shedding 
=  1,  nonsymmetric  shedding 

Separation  flag. 

=  0,  no  separation 
=  1,  laminar  separation 
=  2,  turbulent  separation 

Reverse  flow  separation  flag. 

=  0,  no  reverse  flow  separation 
( preferred ) 

=  1,  reverse  flow  separation 

Vortex  core  model  flag. 

=  0,  potential  vortex  model 
=  1,  diffuse  vortex  model  (pre¬ 
ferred  ) 

Two-dimensional  unsteady  pressure 
term  flag. 

=  0,  omit  3$/3t  term  in  Cp  equation 
=  1,  include  3$/3t  term  (preferred) 

Nose  force  flag. 

=  0,  estimate  nose  forces  ahead  of 
the  first  axial  station  (pre¬ 
ferred  ) 

=  1 ,  zero  nose  forces  ahead  of  the 
first  axial  station 
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NXFV  Number  of  x  stations  at  which  the 

flow  field  is  calculated  and/or  spe¬ 
cial  output  is  generated. 

(0  <:  NXFV  <  8) 

NFV  Number  of  field  points  for  flow 

field  calculation. 

(0  <  NFV  <  200) 

NVP  Number  of  +r  vortices  on  +y  side  of 

body  for  a  restart  calculation. 

(0  <  NVP  <  70) 

NVR  Number  of  reverse  flow  vortices  on 

+y  side  of  body  for  a  restart  calcu¬ 
lation.  (0  <  NVR  <  30) 

NVM  Number  of  -r  vortices  on  -y  side  of 

body  for  a  restart  calculation. 

(0  <  NVM  <  70) 

NVA  Number  of  reverse  flow  vortices  on 

-y  side  of  body  for  a  restart  calcu¬ 
lation.  (0  <  NVA  <  30) 

NASYM  Asymmetric  vortex  shedding  flag. 

=0,  no  asymmetry 
=  1,  forced  asymmetry 

Item  2  is  a  single  card  containing  10  integer  flags/parame¬ 
ters.  These  values  are  used  to  specify  several  options  available 
in  VTXCLD. 
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Item 


Variable 


Description 


2  (1015) 


NHEAD  Number  of  title  cards  in  Item  3. 

( NHEAD  >  0) 

NPRNTP  Pressure  distribution  print  flag. 

=  0,  pressure  distribution  output 
only  at  special  output  sta¬ 
tions  ( NXFV  >  0) 

=  1,  pressure  distribution  output 
at  each  x  station 


NPRNTS  Vortex  separation  print  flag. 

=0,  no  separation  output 

=  1,  separation  point  summary  (pre 
ferred) 

=  2,  detailed  separation  calcula¬ 
tion  (for  diagnostic  purposes 
only) 

NPRNTV  Vortex  field  summary  output  flag. 

=  0,  vortex  field  is  output  at 

special  output  stations  only 
(NXFV  >  0) 

=  1,  vortex  field  is  output  at 
every  x  station 

NPLOTV  Vortex  field  printer-plot  flag. 

=  0,  no  plots 

=  1,  plot  full  cross  section,  con¬ 
stant  scale 

=  2,  plot  half  cross  section,  con¬ 
stant  scale 
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NPLOTA 


NTH 


NCORE 


NSKIN 


NCOMP 


=  3,  plot  full  cross  section,  vari¬ 
able  scale 

Additional  output  flag. 

=  0,  no  additional  output 
=  1,  additional  output  only  at 

special  x  station  (NXFV  >  0) 

=  2,  additional  output  at  each 
station 

Numerical  conformal  mapping  angle 
f  lag . 

=  0,  calculate  circle  0's  for  map¬ 
ping  at  5°  increments  (pre¬ 
ferred  ) 

=  1,  input  6's 

Vortex  core  radius  flag. 

=  0,  maximum  allowable  local  vortex 
core  size  is  . 05*SD 
=  1,  no  upper  limit  on  the  local 

vortex  core  size,  the  vortex 
core  size  is  given  by:  RCORE*D 

Axial  skin  friction  estimation  flag. 
=  0,  no  skin  friction  estimation 
=  1,  axial  skin  friction  estimated 

Compressibility  correction  flag. 

=  0,  incompressible  calculation 
=  1,  compressibility  correction 

applied  (Mach  number  must  be 
input  in  Item  6) 
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Item  3  is  a  set  of  NHEAD  title/summary  cards  which  are 
printed  at  the  top  of  the  first  output  page.  They  are  reproduced 
in  the  output  just  as  they  are  input. 


Item 

Variable 

Description 

3  (20A4) 

TITLES 

NHEAD  title/summary  cards 
output  on  first  page. 

to  be 

Item  4  is 

needed  to  form 

a  single  card  containing  reference  information 
the  aerodynamic  coefficients. 

Item 

Variable 

Description 

4  (5F10.5) 

REFS 

Reference  area,  Sref.  REFS  > 

0. 

REFL 

Q 

Reference  length,  ref.  REFL 

>  0. 

XM 

Axial  position  of  the  moment 

measured  from  the  body  nose. 

center. 

SL 

Body  length,  L. 

SD 

Maximum  body  diameter,  D. 

Item  5  is  a  single  card  containing  the  flow  conditions. 
Item  Variable  Description 

5  (3F10.5)  ALPHAC  Angle  of  incidence  (degrees). 

(0°  <  a  <  90°) 
c 

PHI  Angle  of  roll,  <j>,  (degrees). 
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RE 


Reynolds  number  based  on  maximum 
body  diameter,  SD,  and  V  .(V  D/v) 

■*  OC  CC 


Omit  Item 

Item 

number  M  . 

or 

nonzero . 

I  tern 

6  (F10.5 

Item 

the  run, 
the  vorte> 

I  tern 

7  (8F10 


6  if  NCOMP  =  0. 

6  is  a  single  card  which  contains  the  freestream  Mach 
This  item  is  included  only  if  NCOMP  in  Item  2  is 


Variable  Description 

)  MACH  Freestream  Mach  number. 

7  is  a  single  card  which  specifies  the  axial  extent  of 
the  transition  region,  and  parameters  associated  with 
wake . 

Variable  Description 


.5)  XI 


Initial  x  station,  XI  >  0,  dimen¬ 

sional  quantity.  (Note:  XI  >  DX/2) 


XF 


Final  x  station,  XF  >  XI,  dimen¬ 
sional  quantity. 


DX  x  increment  for  vortex  shedding 

calculation.  Typical  value, 

DX  =  D/2. 


XTRl  Beginning  of  transition  region, 

dimensional  quantity. 

XTR2  End  of  transition  region,  dimen¬ 

sional  quantity. 
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EMKF  Minimum  distance  of  shed  vortex 

starting  position  from  the  body 
surface . 

=  1.0,  shed  vortices  located  ac¬ 
cording  to  stagnation  point 
criterion 

>  1.0,  minimum  radii  position  of 

shed  vortices.  Typical 
value,  EMKF  =  1.05 

RGAM  Vortex  combination  factor. 

=  0.0,  vortices  are  not  combined 
( preferred ) 

>  0.0,  radial  distance  within  which 

vortices  are  combined 


VRF  Vortex  reduction  factor,  6^  to 

account  for  observed  decreases  in 
vortex  strength. 

=  0.6,  for  bodies  with  bases 
=  1.0,  for  closed  bodies 

Item  8  is  a  single  card  containing  the  integration  toler¬ 
ance,  asymmetry  perturbation  information,  alternate  separation 
criteria,  and  the  core  radius. 


I  tern 


Variable 


Description 


8  (8F10.5)  E5  Error  tolerance  for  the  vortex  tra¬ 

jectory  calculation.  Typical  range, 
E5  =  0.01  to  0.05  or  larger  if  ap¬ 
propriate  . 
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XTABL 


Location  at  which  a  table  of  points 
on  the  body  and  on  rhe  circle  are 
printed . 

=  0.0,  no  table  of  points  is 
printed  (preferred) 

>  0.0,  table  of  points  for  all 
x  <  XTABL 

XASYMI  Initial  x  station  for  asymmetric 

perturbation . 

XASYMF  Final  x  station  for  asymmetric  per¬ 

turbation  . 

DBETA  Anqular  displacement  of  separation 

points  for  asymmetry  perturbation. 

SEPL  Stratford  laminar  criterion. 

=  0.0,  program  uses  SEPL  =  .087 
( preferred ) 

>  0.0,  program  uses  input  value 

SEPT  Stratford  turbulent  criterion. 

=  0.0,  program  uses  SEPT  =  .350 
( preferred ) 

>  0.0,  program  uses  input  value 

RCORE  Ratio  of  the  local  core  radius  to 

the  local  body  diameter.  Default 
value  is  .025.  Maximum  allowable 
value  is  .025*D/rQ. 

Item  9  is  a  single  card  containing  two  flags  for  the  source 
distribution.  The  source  distribution  may  be  input  ( NSOR  >  0), 
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calculated  from  the  body  geometry  (NSOR  =  -1),  or  emitted  from 
the  solution  (NSOR  =  0) .  If  the  Mach  number  is  not  equal  to  zero 
(NCOMP  >  0)  then  a  different  source  distribution  is  required  for 
each  Mach  number. 

Item  Variable  Description 

9  (215)  NSOR  Source  distribution  flag. 

=  -1 ,  calculate  the  source  distri¬ 
bution  from  the  input  geome¬ 
try  input 

=0,  no  sources  (not  recommended) 

>  0,  number  of  sources  to  be 

input 

NPRT  Source  distribution  print  flag. 

=  0,  output  source  locations  and 

strengths 

=  1,  above  output  plus  input  and 

source  geometry  comparison 

=  2,  above  output  plus  source- 

induced  surface  velocities 

Omit  Items  10  and  11  if  NSOR  <  0. 

Items  10  and  11  contain  the  nondimensional  source  locations 
and  nondimensional  source  strengths,  respectively. 

I tern  Variable  Description 


10  (6E12.5)  XSRC  Nondimensional  source  locations, 

x/L.  (NSOR  values,  6  per  card) 
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Item 


Variable 


Description 


11  (6E12.5)  QS  Nondimensional  source  strengths, 

2 

Q/4ttL  V  cosa  . 

op  c 

Items  12  through  15  contain  the  body  geometry  information. 
This  information  must  be  included  for  all  bodies  whether  in  cir¬ 
cular  or  noncircular. 

Item  Variable  Description 

12  ( I 5 )  NXR  Number  of  entries  in  the  body  geome¬ 

try  table.  (1  <  NXR  <  101) 


Item 


Variable 


Description 


13  (8F10.5)  XR  Nondimensional  axial  stations  in  the 

geometry  table,  x/L.  (NXR  values,  8 
per  card ) 


Item 


Variable 


Description 


14  (8F10.5)  R  Nondimensional  body  radius  at  the 

NXR  x/L-stations,  rQ/L  (NXR  values, 
8  per  card)  .  For  an  ellipse,  rQ  = 
(a  +  b)/2.  For  arbitrary  bodies,  rQ 
is  equivalent  radius. 
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Item 


Variable 


Items  16  and  17  contain  the  nondimensional  horizontal  and 
vertical  half  axis  lengths  if  the  cross  sections  are  elliptical, 
NCIR  =  1. 


Item 


Variable 


Description 


16  (8F10.5)  AE  Nondimensional  horizontal  half  axis 

length  of  the  elliptical  cross  sec¬ 
tion  at  the  NXR  x/L-stations,  a/L. 
(NXR  values,  8  per  card) 


Item 


Variable 


Description 


17  (8F10.5)  BE  Nondimensional  vertical  half  axis 

length  of  the  elliptical  cross  sec¬ 
tion  at  the  NXR  x/L-stations,  b/L. 
(NXR  values,  8  per  card) 


Omit  Items  18  through  29  if  NCIR  2. 

Items  18  through  29  are  used  to  specify  a  body  of  arbitrary 
cross  section  shape.  This  body  should  also  be  polygonal  in  shape 
(i.e.,  negative  interior  angle  can  cause  problems  in  the  numer¬ 
ical  conformal  mapping  procedure). 
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Omit  Item  18  if  NTH  =  0 


Item 

18  (8F10 


Item 

19  (215) 


I  tern 

20  (8F10 

Omit  Items 


Variable  Description 

5)  TH  Special  table  of  circle  angles  for 

arbitrary  body  mapping,  (i  =  1,  73), 
or  to  be  used  only  if  the  density  of 
points  on  the  cross  section  needs 
adjusting  to  improve  the  resolution 
of  pressure  distribution.  Not  rec¬ 
ommended  under  normal  conditions. 


Variable  Description 

MNFC  Number  of  mapping  coefficients  for 

the  arbitrary  body.  (MNFC  100) 

MXFC  Number  of  axial  stations  at  which 

the  arbitrary  body  is  specified. 
(1  <  MXFC  <  10).  For  a  similar 
shape  at  all  cross  sections,  MXFC=1. 


Variable  Description 

.5)  XFC  x  stations  at  which  the  arbitrary 

body  is  specified.  (MXFC  values,  8 
per  card).  For  a  similar  cross 
section  body,  MXFC=1,  XFC(l)  <  XI. 

21  through  23  if  NCF  >  0. 
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Items  21  through  23  are  repeated  for  each  of  the  MXFC  axial 
stations  when  the  transformation  coefficients  are  to  be  calcu¬ 
lated  (NFC  =0). 


Item 


Variable 


Description 


21  (15)  NR  Number  of  points  describing  the 

arbitrary  body  cross  section 
(2  <  NR  <  30).  (Note:  NR  must  be 
the  same  for  all  x  stations. ) 


22  (6E12.5)  XRC  Horizontal  coordinates  of  the  arbi¬ 

trary  body  cross  section.  (NR 
values,  6  per  card).  The  convention 
for  ordering  the  coordinates  from 
bottom  to  top  in  a  counter-clockwise 
fashion  is  observed,  Figure  9. 
Right/left  body  symmetry  is  re- 
guired . 


23  (6E12.6)  YRC  Vertical  coordinates  of  the  arbi¬ 

trary  body  cross  section.  (NR 
values,  6  per  card).  See  Figure  9. 


Omit  Items  24  through  29  if  NCF  =  0. 


Items  24  through  29  are  repeated  for  each  of  the  MXFC  axial 
stations  if  the  numerical  transformations  are  available  from  a 
previous  calculation.  The  only  purpose  for  providing  the  option 
to  input  the  mapping  parameters  is  to  eliminate  the  need  to  re- 
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compute  the  numerical  mapping  in  subsequent  runs  and  save  compu¬ 
tation  time. 

Item  Variable  Description 

24  (2F10.5)  ZZC  Mapping  offset  distance. 

RZC  Equivalent  radius  of  transformed 

cross  section. 


Item 


Variable 


Description 


25  (15)  NR  Number  of  points  describing  the 

arbitrary  body  cross  section 
(NR  <  30).  (Note:  NR  must  be  the 
same  for  all  x  stations.) 


Item 


Variable 


Description 


26  (6E12.5)  XRC  Horizontal  coordinates  of  the  arbi¬ 

trary  body  cross  section.  (NR 
values,  6  per  card).  See  Figure  9 
and  convention  for  Item  22. 


Item 


Variable 


Description 


27  (6E12.5)  YRC  Vertical  coordinates  of  the  arbi¬ 

trary  cross  section.  (NR  values,  6 
per  card).  See  Figure  9. 
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Item 


Variable 


Description 


28  (6E12.5)  THC  Circle  angles  6,  of  body  points. 

(NR  values,  6  per  card) 

Item  Variable  Description 

29  (6E12.5)  AFC  Mapping  coefficients.  (MNFC  values, 

6  per  card) 

Omit  Items  30  and  31  if  NXFV  =  0. 

Item  30  is  a  single  card  containing  the  x  stations  at  which 
additional  output  is  printed  and/or  calculated  if  NXFV  >  0. 


Item  Variable  Description 

30  (8F10.5)  XFV  x  stations  at  which  additional  out¬ 

put  is  printed  or  calculated.  (NXFV 
values,  8  values  maximum) 

Omit  Item  31  if  NFV  =  0. 

Item  31  is  a  set  of  NFV  cards  which  contains  the  coordinates 
of  the  field  points  where  the  velocity  components  are  to  be  cal¬ 
culated  at  each  of  the  axial  stations  in  Item  30. 


Item  Variable  Description 

31  (2F10.5)  YFV,ZFV  y  and  z  coordinates  of  the  field 

points  for  the  velocity  field  calcu¬ 
lation.  (NFV  cards,  maximum  of  200) 
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It  is  important  that  all  points  lie 
outside  the  body  surface  at  all 
axial  stations. 

Omit  Items  32  through  37  if  NVP  +  NVR  +  NVM  +  NVA  =  0 

Items  32  through  36  provide  information  needed  to  restart 
the  calculation  procedure.  It  is  not  necessary  to  input  all 
four  types  of  vorticity  in  a  restart  calculation. 

Item  32  is  a  single  card  containing  the  nose  force  and  mo¬ 
ment  coef t icients  at  the  restart  point;  that  is,  X  =  XI .  These 
values  may  be  set  equal  to  zero,  but  the  resulting  forces  and 
moments  from  the  current  calculation  will  be  in  error. 


Item  Variable 


32  (5F10.5)  CN 

CY 

CA 

ODI 


CM 

CR 

CSL 


Description 

Normal-force  coefficient 

Side-force  coefficient 

Axial  pressure  force  coefficient 

Axial  skin  friction  force  coeffi¬ 
cient 

Pitching-moment  coefficient 
Yawing-moment  coefficient 
Rolling -moment  coefficient 
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Item  33  is  a  set  of  NVP  cards  which  specify  a  cloud  of  posi¬ 
tive  primary  body  separation  vorticity,  usually  on  the  right  side 
of  the  body.  The  variable  XSHEDP  associated  with  each  vortex  is 
used  only  to  identify  individual  vortices  and  permit  the  user  to 
follow  individual  vortices  during  the  calculation. 


Omit  Items  33  if  NVP 

I  tern  Variable 

33  (4F10.5)  GAMP 

YGP 

ZGP 

XSHEDP 


0. 

Description 

+y  primary  vorticity,  r/V^. 
Horizontal  coordinate,  y,  of  vortex. 
Vertical  coordinate,  z,  of  vortex. 
Origin  of  the  vortex. 


Item  34  is  a  set  of  NVR  cards  which  specify  a  cloud  of  re¬ 
verse  flow  (secondary)  body  separation  vorticity.  This  array  is 
a  convenient  way  to  put  an  arbitrary  cloud  of  additional  vorti¬ 
city  which  is  to  be  maintained  separately  from  the  other  body 
vorticity  in  the  field. 


Omit  Item  34  if  NVR  =  0. 


I  tern 

Variable 

Description 

34  (4F10.5) 

GAMR 

+y  reverse  flow  vorticity. 

r/v_. 

YGR 

Horizontal  coordinate,  y, 

of  vortex 

ZGR 

Vertical  coordinate,  z,  of 

vortex . 
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XSHEDR 


Origin  of  the  vortex 


Item  35  is  a  set  of  NVM  cards  which  specify  a  cloud  of  nega¬ 
tive  primary  body  separation  vorticity,  usually  on  the  left  side 
of  the  body.  In  the  case  of  a  symmetric  flow  field  (ISYM^O), 
this  block  of  vorticity  is  automatically  defined  as  the  mirror 
imaqe  of  the  positive  vorticity  input  in  Item  33. 

Omit  Item  35  if  ISYM  =  0  or  NVM  =  0. 


I  tem 

Variable 

Description 

35  (4F10.5) 

GAMM 

-y  primary  vorticity,  r/V^. 

YGM 

Horizontal  coordinate,  y,  of  vortex. 

ZGR 

Vertical  coordinate,  z,  of  vortex. 

XSHEDR 

Origin  of  the  vortex. 

Item  36  is  a  set  of  NVA  cards  which  specify  a  cloud  of  re¬ 
verse  flow  (secondary)  body  separation  vorticity.  This  block  of 
vorticity  is  analogous  to  Item  34,  and  it  is  omitted  if  NVA  =  0 
or  if  ISYM  =  0.  As  with  the  previous  item,  this  block  of  vorti¬ 
city  is  automatically  defined  as  the  mirror  image  of  the  negative 
vorticity  input  as  Item  34  if  symmetry  is  required  by  ISYM  =  0. 

Omit  Item  36  if  ISYM  =  0  or  NVA  =  0. 


Item 

Variable 

Description 

36  (4F10.5) 

GAMA 

-y  reverse  flow  vorticity,  r/V  . 

-4  OP 

YGA 

Horizontal  coordinate,  y,  of  vortex 
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ZGA 


Vertical  coordinate,  z,  of  vortex. 


ZSHEDA  Origin  of  the  vortex. 

Item  37  contains  the  pressure  distribution  at  the  initial  x 
station  for  a  restart  calculation.  If  a  symmetric  case  is  being 
considered,  ISYM=0,  only  half  of  the  pressure  distribution  is 
needed. 


Item 


Variable 


37  (3F10.5)  YBO 

ZBO 

CPO 


Description 

y  body  coordinate. 

z  body  coordinate. 

pressure  coefficients  at 
(YBO,  ZBO). 


This  concludes  the  input  deck  for  a  single  run  of  VTXCLD. 
The  code  is  not  set  up  to  stack  input  for  multiple  cases  because 
of  the  sometimes  long  and  generally  unpredictable  run  times. 

Input  Preparation 

In  this  section,  some  additional  information  is  provided  to 
assist  the  user  in  the  preparation  of  input  for  certain  selected 
problems.  The  previous  section  on  the  input  description  must  be 
used  to  understand  the  individual  variables  which  go  into  VTXCLD, 
and  this  section  will  permit  the  user  to  select  the  appropriate 
input  to  get  optimum  results  from  the  code. 

Numerical  mapping.-  The  specification  of  the  appropriate 
numerical  mapping  parameters  (Items  18  through  29)  depends  on  the 
shape  of  the  cross  section.  In  the  interest  of  optimum  computa- 
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tion  time,  the  fewest  possible  terms  in  the  transformation  series 
( MNFC  in  Item  19)  should  be  used  which  will  produce  an  adequate 
mapping.  Each  different  shape  should  be  checked  by  the  user  to 
determine  the  proper  number  of  coefficients.  For  example,  a 
simple  cross  section  which  is  similar  to  a  circle  may  require  as 
few  as  10  coefficients.  Others  may  require  as  many  as  20  coeffi¬ 
cients  for  an  excellent  mapping.  If  a  large  number  of  calcula¬ 
tions  are  to  be  made  for  a  noncircular  shape,  it  is  worthwhile 
for  the  user  to  try  different  numbers  of  coefficients  in  an  ef¬ 
fort  to  find  an  optimum  number.  Not  only  should  the  shapes  be 
compared  carefully,  the  unseparated  flow  pressure  coefficients 
should  be  compared  to  evaluate  the  quality  of  the  mapping. 

The  numerical  mapping  now  included  in  VTXCLD  has  been  used 
for  a  variety  of  cross  section  shapes,  and  very  few  problems  have 
been  observed.  One  general  problem  area  is  bodies  with  concave 
sides.  The  mapping  procedure  has  difficulty  converging  on  a  set 
of  coefficients.  Usually  the  body  shape  can  be  modified 
slightly,  to  give  the  concave  region  a  flat  or  slightly  convex 
shape,  and  the  mapping  converges  easily.  This  small  change  in 
the  shape  probably  has  very  little  effect  on  the  predicted  pres¬ 
sure  distribution.  This  is  an  example  of  the  type  of  problem  to 
which  the  user  should  be  alert. 

Integration  interval.-  Choice  of  the  appropriate  integration 
interval  ( DX  in  Item  7)  depends  to  a  large  extent  on  the  type  of 
results  of  interest.  If  the  user  is  interested  in  detailed  nose 
loads  and  separation  on  the  nose,  where  the  radius  is  changing 
rapidly,  the  internal  should  be  smaller  than  the  recommended 
D/2.  On  very  long  bodies,  the  recommended  interval  is  usually 
adequate  for  accurate  predictions. 

Vortex  core.-  The  core  size  specification  is  referenced  to 
the  local  body  diameter  (RCORE  in  Item  8).  The  default  value  is 
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0.025,  and  this  value  is  adequate  for  most  cases.  The  larger  the 
core  radius  specified,  the  more  the  vortex  effects  are  reduced. 
To  eliminate  the  possibility  of  negating  the  vortex  effect  en¬ 
tirely,  a  maximum  vortex  core  of  .05  times  the  maximum  body  di¬ 
ameter  is  built  into  the  code.  Except  in  unusual  situations,  the 
default  value  is  recommended,  and  this  is  the  value  used  in  the 
comparisons  with  data  shown  in  a  later  section. 

Sample  Cases 

In  this  section,  nine  sample  input  cases  are  described  to 
illustrate  the  various  program  options  available  in  VTXCLD.  The 
purpose  is  to  point  out  the  available  features  and  options  and  to 
provide  a  range  of  sample  inputs  which  will  help  users  prepare 
input  for  specific  cases.  Table  I  provides  a  summary  of  the 
sample  input  cases  to  assist  the  user  in  finding  an  appropriate 
example . 

Sample  Case  1,  Figure  10(a),  is  a  test  run  based  on  the  3- 
caliber  ogive-cylinder  body  in  Reference  20.  This  run  will  also 
be  used  as  the  sample  output  case  in  the  next  section.  It  is 
recommended  that  this  case  be  run  initially  to  provide  a  check 
with  the  results  presented  in  a  following  section.  The  input  for 
the  first  sample  case  is  described  in  detail,  and  for  the  fol¬ 
lowing  cases,  only  the  major  differences  are  discussed. 

Input  Item  1  indicates  that  the  case  consists  of  a  body  with 
a  circular  cross  section,  NCIR  =  0,  and  a  symmetric  vortex  wake, 
ISYM  =  0.  The  separation  criteria  is  laminar,  NBLSEP  =  1,  and  a 
diffuse  vortex  core  model  is  to  be  used,  NDFUS  =  1.  NDPHI  =  1 
indicates  that  the  two-dimensional  unsteady  pressure  term  is 
included  in  pressure  calculation.  NXFV  =  7  and  NFV  =  0  are  in¬ 
cluded  to  denote  special  output  to  be  printed  for  seven  axial 
stations  which  are  specified  in  Item  30.  If  NFV  >  0,  the  flow 
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field  will  be  calculated  at  the  7  axial  stations  in  addition  to 
any  special  output. 

Input  Item  2  indicates  there  will  be  four  title  cards  input 
in  Item  3.  NPRNTP  =  0  specifies  the  pressure  distribution  be 
output  at  the  seven  special  x  stations  specified  by  NXFV,  and 
NPRNTS  =  1  causes  a  summary  of  the  separation  points  to  be  output 
at  all  x  stations.  The  flag  NPRNTV  =  0  results  in  a  summary  of 
the  vortex  field  at  the  seven  axial  stations  specified.  NPLOTV  = 
3  and  NPLOTA  =  1  specifies  that  a  line  printer  plot  of  the  vortex 
field  is  created  at  the  seven  special  x  stations  if  separation 
has  occurred.  NSKIN  =  1  will  result  in  an  estimation  of  the  skin 
friction  on  the  body,  and  NCOMP  =  0  specifies  an  incompressible 
calculation . 

Input  Item  3  is  a  set  of  four  title/summary  cards  which  are 
printed  at  the  top  of  the  first  output  page. 

Input  Item  4  contains  the  reference  information  for  the 
run.  The  reference  area,  REFS,  is  the  maximum  cross  sectional 
area,  the  reference  length,  REFL ,  is  the  maximum  diameter,  and 
the  moment  center,  XM,  is  located  at  the  end  of  the  ogive  nose. 
The  length  of  the  body,  SL,  is  50.478  inches,  and  the  maximum 
body  diameter,  SD,  is  4.7  inches. 

Item  5  contains  the  flow  conditions.  The  angle  of  attack, 
ALPHAC  is  15  degrees  and  the  Reynolds  number,  RE,  is  440000. 
There  is  no  roll  angle  for  this  case,  PHI  =  0.0. 

Item  7  specifies  run  parameters.  The  calculation  procedure 
begins  at  x  station,  XI,  2.35  and  continues  to  station,  XF,  47.0 
in  DX  increments  of  2.35  or  one  body  radius.  Separation  vortices 
are  initially  placed  a  minimum  1.05  radii,  EMKF,  from  the  axis  of 
the  missile.  The  vortex  reduction  factor,  VRF ,  is  set  to  0.6  as 
recommended  for  bodies  which  do  not  close  themselves. 
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The  only  major  item  of  interest  in  input  Item  8  is  the  tol¬ 
erance  for  the  vortex  trajectory  calculation,  E5.  This  is  set  to 
0.05,  a  typical  value.  Some  additional  comments  on  this  tol¬ 
erance  are  contained  in  the  previous  section  on  Error  Messages 
and  Stops. 

Input  Item  9  contains  the  flags  for  the  source  distribution 
calculation.  NSOR  =  -1  specifies  the  source  distribution  is 
unknown  and  is  to  be  calculated  by  VTXCLD.  NPRT  =  2  provides  for 
additional  output  regarding  the  source  distribution.  This  output 
includes  a  comparison  of  the  input  geometry  and  the  geometry 
calculated  from  the  source  distribution  which  provides  a  good 
means  of  checking  the  calculated  source  distribution.  Also 
printed  are  the  surface  velocities  calculated  by  the  source  dis¬ 
tribution. 

Input  Item  12  contains  NXR,  the  number  of  entries  in  the 
body  geometry  table.  Items  13,  14,  and  15  are  the  axial  sta¬ 
tions,  XR,  the  corresponding  body  radii,  R,  and  the  body  slopes, 
DR,  respectively.  The  axial  stations  and  the  radii  are  nondimen- 
sionalized  with  respect  to  the  body  length,  SL. 

The  final  input  for  this  case  is  Item  30  which  identifies 
XFV,  the  seven  axial  stations  where  special  output  is  printed. 

The  following  sample  case  inputs  are  discussed  only  briefly 
with  some  elaboration  of  the  parameters  of  special  interest. 

Sample  Case  2,  Figure  10(b),  is  the  same  as  Sample  Case  1 
except  the  source  distribution  is  input  and  the  angle  of  attack 
is  changed.  NSOR  =  61  in  Item  9  indicates  that  there  are  61 
source  locations,  XSRC,  and  strengths,  OS,  input  via  Items  10  and 
11,  respectively.  The  locations  and  strengths  were  obtained  from 
the  results  of  Sample  Case  1.  The  angle  of  attack,  ALPHAC,  in 
Item  5  is  changed  to  20  degrees. 
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Sample  Case  3,  Figure  10(c),  is  a  body  with  2:1  elliptical 
cross  section.  This  model  is  described  in  Reference  21.  The 
elliptical  cross  section  is  handled  by  an  analytical  conformal 
mapping,  NCIR  =  1  in  Item  1,  and  the  other  inputs  required  for 
the  ellipse  are  the  half-axis  lengths,  AE  and  BE,  input  in  Items 
16  and  17. 

Sample  Case  4,  Figure  10(d),  is  a  body  of  square  cross  sec¬ 
tion  with  rounded  corners  from  Reference  22.  A  compressible  Mach 
number  is  specified.  For  this  cross  section,  a  numerical  con¬ 
formal  mapping  must  be  used,  NCIR  =  2  in  Item  1.  NCF  =  0  in  Item 
1  indicates  that  the  mapping  coefficients  are  to  be  calculated. 
The  desired  number  of  mapping  coefficients,  MNFC,  is  set  to  20  in 
Items  19.  Item  19  also  specifies  the  number  of  axial  stations  at 
which  the  body  is  defined,  MXFC  =  1.  MXFC  =  1  is  used  because 
all  of  the  axial  stations  are  similar  in  shape;  therefore,  it  is 
necessary  to  calculate  mapping  coefficients  only  at  a  single  sta¬ 
tion.  If  the  cross  sections  along  the  body  are  not  similar, 
mapping  coefficients  must  be  calculated  at  a  number  of  different 
sections.  Item  20  is  the  axial  station,  XFC,  at  which  the  map¬ 
ping  coefficients  are  calculated.  In  Item  21,  NR  =  25  is  the 
number  of  points  describing  the  body  cross  section  and  Items  22 
and  23  contain  the  horizontal  and  vertical  coordinates  of  the 
cross  section,  respectively. 

Compressibility  is  included  by  setting  the  compressibility 
flag,  NCOMP,  to  1  in  Item  1.  The  freestream  Mach  number  is 
defined  in  Item  6,  MACH  =  0.5. 

Sample  Case  5,  Figure  10(e),  illustrates  the  required  input 
for  an  axisymmetric  body  at  nonzero  Mach  number.  This  case  is 
Sample  Case  1  with  a  compressibility  correction.  NCOMP  in  Item  2 
is  set  to  1,  and  the  Mach  number,  MACH,  is  set  to  0.8  in  Item 
6.  The  source  distribution  is  calculated  because  it  changes  with 
Mach  number. 
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Sample  Case  6,  Figure  10(f),  illustrates  the  specification 
of  a  transition  region.  The  beginning  and  end  of  the  transition 
region  are  specified  in  Item  7.  The  beginning,  XTR1  =  7.05,  and 
the  end,  XTR2  =  11.75,  define  the  region  over  which  the  separa¬ 
tion  criteria  changes  from  laminar  to  turbulent. 

Sample  Case  7,  Figure  10(g),  illustrates  the  reguired  input 
for  a  case  in  which  an  asymmetric  wake  occurs.  In  order  to  ob¬ 
tain  an  asymmetric  vortex  wake,  the  flags  ISYM  and  NASYM  in  Item 
1  are  set  to  1.  The  reguired  symmetric  perturbation  is  specified 
in  Item  8  where  the  initial  x  station  for  the  perturbation, 
XASYMI  =  0.0,  and  the  final  station,  XASYMF  =  5.145,  are  de¬ 
fined.  The  angular  displacement  of  the  separation  points,  DBETA , 
is  2.0  degrees  over  the  entire  region.  In  Item  2,  NPRNTP  and 
NPRNTV  have  been  set  to  1.  This  will  cause  the  pressure  distri¬ 
bution  and  vortex  wake  to  be  printed  at  each  x  station. 

The  velocity  field  is  calculated  at  57  (NFV)  field  points  at 
seven  axial  stations,  NXFV.  These  inputs  are  defined  in  Item 
1.  The  seven  axial  stations  where  the  velocity  field  is  desired 
are  specified  in  Item  30,  and  the  57  field  points  where  the  ve¬ 
locities  are  calculated  are  included  in  Item  31. 

Sample  Case  8,  Figure  10(h),  is  a  lobe  body  from  Reference 
23  shown  to  illustrate  another  arbitrary  cross  section  body.  The 
noncircular  input  has  similar  format  to  the  rounded  sguare  body 
of  Sample  Case  4,  and  the  only  difference  is  in  the  body  cross 
section  specification  in  Items  22  and  23.  No  source  distribution 
is  input  or  calculated  for  this  case;  therefore,  this  situation 
approximates  a  two-dimensional  example  in  which  the  axial  coordi¬ 
nate  represents  time. 

Sample  Case  9,  Figure  10(i),  is  an  example  of  a  restart 
calculation  for  Sample  Case  2.  The  case  is  restarted  from  XI  = 
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21.15  and  continued  to  XF  =  47.0.  Additional  input  required  for 
this  case  are  the  overall  loadings  in  Item  32,  the  vortex  loca¬ 
tions  and  strengths  in  Items  33  and  35,  and  the  pressure  distri¬ 
bution  in  Item  37.  This  case  is  included  to  illustrate  the  means 
to  restart  a  calculation  which  did  not  finish  or  one  in  which  the 
vortex  field  was  adjusted  to  include  additional  vortices  from  a 
forward  lifting  surface. 

Output  Description 

A  typical  output  file  from  VTXCLD  is  described  in  this  sec¬ 
tion.  In  general,  the  output  quantities  from  VTXCLD  are  labeled 
and  each  page  is  headed  with  approximate  descriptive  informa¬ 
tion.  The  output  from  Sample  Case  1,  Figure  10(a),  requires  52 
pages,  and  it  will  be  described  briefly  page  by  page.  Represen¬ 
tative  pages  are  shown  in  Figure  11. 

The  first  output  page  contains  the  title/summary  cards,  the 
reference  quantities,  the  flow  conditions,  the  initial  run  condi¬ 
tions,  and  the  run  flags.  Page  1  is  simply  an  echo  of  the  input 
parameters.  Also  contained  on  this  page  are  the  force  and  moment 
definitions.  Pages  2  and  3  echo  the  input  body  geometry,  XR,  R, 
and  DR.  The  dimensional  X  and  R  values  are  also  printed. 

Pages  4  through  10  contain  information  from  the  three-dimen¬ 
sional  source  distribution  calculation.  Page  4  is  a  summary  of 
the  source  placement  and  the  body  geometry  at  these  locations. 
Page  5  shows  the  source  locations  and  strengths,  and  pages  6  and 
7  are  a  comparison  of  the  input  geometry  shape,  R/L(INPT),  and 
the  shape  calculated  from  the  source  distribution  R/L(S.D.). 
This  comparison  is  indicative  of  the  quality  of  the  source  dis¬ 
tribution.  A  second  indication  of  difficulties  with  the  source 
distribution  will  show  up  in  the  number  of  iterations  required  to 
find  the  proper  shape,  M-INDEX.  If  this  number  is  less  than  four 
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or  five,  the  chances  are  good  the  source  distribution  is  a  good 
representation  of  the  body.  If  this  number  is  large,  examine  the 
input  and  predicted  shapes  carefully  for  problems.  Pages  8  and  9 
are  the  surface  velocities  calculated  by  the  source  distribution, 
and  they  should  be  examined  to  ensure  that  no  exceedingly  large 
velocities  occur,  another  indication  of  a  problem  with  the  body 
model.  Small  oscillations  in  the  velocities  on  the  cylindrical 
portion  of  a  body  are  to  be  expected. 

Page  10  contains  a  summary  of  source  locations  and 
strengths,  and  the  values  on  this  page  are  in  the  correct  format 
for  input  into  VTXCLD  for  additional  cases  using  the  same  geome¬ 
try.  See  Sample  Case  2  for  an  example.  Also  contained  on  page 
10  are  the  x  stations  at  which  special  output  is  printed. 

Page  11  starts  the  output  from  the  analysis  part  of  the 
code.  The  first  numeric  line  on  page  11  contains  the  axial  loca¬ 
tions,  the  local  radius,  the  body  slope,  and  the  transition 
Reynolds  number,  RETR.  Since  axial  station  x  *  2.35  is  a  special 
output  station  specified  in  Item  30,  the  pressure  distribution  is 
printed.  The  predicted  results  on  pages  11  and  12  are  a  summary 
of  the  velocities  and  pressures  on  the  circumference  at  this 
particular  x  station.  The  coordinates  are  dimensional,  and  the 
velocities  are  made  dimensionless  by  the  freestream  velocity. 
The  pressure  coefficient  is  followed  by  the  unsteady  term,  dij>/dt. 

Following  the  pressure  distribution  is  a  summary  of  the 
section  loadings  and  the  overall  loadings  up  to  the  current  sec¬ 
tion.  CN(X),  CY(X),  and  CA(X)  are  the  normal-,  side-,  and  axial- 
force  coefficients  per  unit  length,  respectively.  CN,  CY,  and  CA 
are  the  overall  normal-,  side-,  and  axial-force  coefficients  from 
the  nose  to  the  current  x  station,  respectively.  CM,  CR,  and  CSL 
are  the  pitching-,  yawing-,  and  rolling-moment  coefficients, 
respectively.  XCPN  and  XCPY  are  the  x  locations  of  the  centers 
of  pressures  for  the  normal  and  side  forces,  respectively. 
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The  next  line  of  output  is  the  estimation  of  the  axial  skin 
friction  coefficient  followed  by  the  total  axial  force  coeffi¬ 
cient  made  up  of  the  pressure  force  and  the  skin  friction  force. 

The  final  block  of  information  on  page  12  contains  a  summary 
of  the  pressure  distribution  and  the  separation  points  at  x  sta¬ 
tion  2.35.  It  is  noted  here  that  no  separation  occurs  at  this 
station. 

Output  pages  13  and  14  have  the  same  format  as  pages  11  and 
12  so  they  will  not  be  discussed  in  detail.  The  x  station  is 
4.7,  a  special  output  station,  so  the  pressure  distribution  is 
again  output.  The  last  two  lines  of  page  14  indicate  that  sepa¬ 
ration  occurs  at  this  axial  station.  The  last  block  of  informa¬ 
tion  contains  the  intitial  locations  and  strengths  of  the  sepa¬ 
rated  vortices. 

Page  15  shows  the  force  and  moment  summary  and  separation 
information  for  x  =  7.05.  This  format  is  the  same  as  the  pre¬ 
vious  stations;  however,  the  pressure  distribution  is  not  output 
at  this  station  because  x  =  7.05  is  not  a  special  output  station. 

Axial  station  9.4  is  again  a  special  output  station.  Page 
16  is  a  summary  of  the  vortex  wake  at  this  station  after  the 
trajectories  of  the  individual  vortices  have  been  calculated. 
Each  individual  vortex  in  the  field  is  described  in  terms  of  its 
location  in  the  real  plane  (Y,  Z)  and  the  transformed  circle 
plane  (YC,  ZC).  The  last  block  of  information  is  the  centroid  of 
vorticity  on  each  side  of  the  body.  This  includes  all  the  vorti- 
city,  including  that  in  the  feeding  sheet;  therefore,  the  numer¬ 
ical  centroid  location  may  be  lower  than  expected  because  of  the 
effect  of  the  vortices  in  the  feeding  sheet. 
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Page  17  has  the  same  format  as  page  11  described  pre¬ 
viously.  Page  18  is  a  line  printer  plot  of  the  cross  section  and 
the  vortex  wake.  Page  19  is  a  summary  of  the  separation  points 
at  x  =  9.4,  as  was  described  for  page  12. 

The  remainder  of  the  output  pages,  except  the  final  three, 
repeat  the  formats  described  previously.  A  typical  output  sta¬ 
tion  has  a  force  and  moment  and  separation  summary  like  page 
15.  Special  output  stations  have  a  wake  summary,  pressure  dis¬ 
tribution  output,  line  printer  plot,  and  separation  summary  like 
pages  16  -  19. 

The  last  three  pages  of  the  output  shown  in  Figure  ll(t), 
(u)  and  (v  give  a  summary  of  the  total  loadings,  the  vortex 
wake,  and  the  pressure  distribution.  Page  50  contains  a  line 
printer  plot  of  the  vortex  wake  at  the  final  x  station  and  a 
summary  of  the  total  loads  on  the  body.  Page  51  is  a  summary  of 
the  vortex  wake  in  tabular  form.  This  information  can  be  used 
for  a  restart  of  the  calculation.  Page  52  is  a  summary  of  the 
pressure  distribution  at  the  last  x  station;  needed  for  a  restart 
calculation . 


RESULTS 


For  purposes  of  evaluating  the  accuracy  and  range  of  appli¬ 
cability  of  the  subsonic  vortex  shedding  model  and  associated 
computer  coce  VTXCLD,  comparisons  of  measured  and  predicted  aero¬ 
dynamic  characteristics  are  presented.  Results  for  both  circular 
and  noncircular  bodies  are  shown  for  a  wide  range  of  flow  condi¬ 
tions  including  flow  asymmetry  and  compressibility  effects. 
Since  the  major  objective  is  the  validation  of  the  prediction 
method,  not  all  predicted  results  are  in  good  agreement  with 
experiment  to  illustrate  to  the  user  where  problems  may  occur 
during  general  use  of  the  code.  Typical  results  from  the  predic¬ 
tion  method  follow. 


-71- 


Circular  Bodies 


The  prediction  method  was  applied  to  an  ogive-cylinder  model 
(Ref.  20)  in  subsonic  flow  over  a  range  of  angles  of  attack  to 
evaluate  the  pressure  prediction  ability  of  the  code.  The  con¬ 
figuration  has  a  three-caliber  ogive  nose  and  a  7.7-caliber  cy¬ 
lindrical  afterbody.  Circumferential  pressure  distributions  at  a 
number  of  axial  stations  are  available  for  a  range  of  flow  condi¬ 
tions. 

As  shown  in  Sample  Case  1,  the  model  is  represented  by  61 
point  sources  and  sinks  on  the  axis,  and  the  ogive  nose  is  used 
to  close  the  model  at  the  base.  Because  of  the  low  Reynolds 
numbers  of  the  experiment,  laminar  separation  is  used  for  all 
predictions  for  this  model.  Comparisons  of  measured  and  pre¬ 
dicted  normal  force  distributions  at  o  =  10,  15,  and  20  degrees 
are  presented  in  Figure  12.  The  predicted  results  with  vortex 
wake  effects  included  are  in  good  agreement  with  experiment  for 
the  angles  shown.  The  significant  influence  of  the  lee  side 
vorticity  is  illustrated  for  a  =  10  and  20  degrees  where  the 
potential  flow  or  slender  body  theory  results  without  vortex 
induced  effects  are  shown  as  a  dashed  curve.  Note  that  the  vor¬ 
tex  effects  are  small  on  the  nose  at  a  =  10  degrees  and  they 
begin  to  grow  as  the  vortex  increases  in  strength  toward  the  end 
of  the  afterbody.  At  a  =  20  degrees,  the  vortex  effects  are 
important  on  the  nose,  and  they  become  a  dominant  effect  on  the 
cylindrical  portion  of  the  model. 

Details  of  the  predicted  characteristics  are  shown  in  Figure 
13  where  comparisons  of  measured  and  predicted  pressure  distribu¬ 
tions  for  a  =  15  and  20  degrees  are  presented.  The  predicted 
results  including  vortex  induced  effects  are  represented  by  the 
solid  curves  and  those  without  vortex  effects  are  shown  as  dashed 
curves.  At  a  =  15  degrees  in  Figure  13(a),  the  vortex  effects  do 
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not  become  significant  until  approximately  4.5  diameters  from  the 
nose.  Some  of  the  roughness  in  the  rational  flow  model  results 
is  caused  by  individual  vortices  moving  too  near  the  body  surface 
during  the  trajectory  calculation.  The  vortex  core  model  does 
tend  to  smooth  the  vortex  induced  effects,  but  there  can  still  be 
a  large  local  effect.  These  local  irregularities  have  a  minimal 
effect  on  the  integrated  loads. 

The  higher  angle  results  shown  in  Figure  13(b)  exhibit  much 
larger  vortex  induced  effects  because  the  separation  vortices  are 
stronger  and  occur  earlier  on  the  body.  The  fact  that  there  are 
more  vortices  with  greater  individual  strength  causes  the  local 
irregularities  to  be  greater.  This  minor  problem  could  be  cor¬ 
rected  to  some  extent  by  decreasing  the  axial  interval  (DX)  of 
separation.  The  individual  strengths  of  the  vortices  would  de¬ 
crease,  but  there  would  be  more  vortices  forming  the  cloud. 
Smoother  pressure  distributions  may  result  at  the  expense  of 
additional  computation  costs. 

The  predicted  vortex  cloud  on  the  ogive  cylinder  model  at 
a  =  15  degrees  is  shown  in  Figure  14.  As  observed  above,  the 
cloud  does  not  become  significant  until  approximately  4.5  diame¬ 
ters  from  the  nose  even  though  separation  begins  at  one  diameter 
from  the  nose. 

The  total  normal-force  coefficient  is  shown  in  Figure  15 
where  both  potential-flow  and  rational-flow  model  results  are 
compared  with  experiment.  As  noted  in  the  pressure  coefficient 
comparisons,  the  vortex  induced  effects  become  significant  at 
a  =  10  degrees  and  grow  with  increasing  angle.  At  low  angles  the 
predicted  results  are  in  good  agreement  with  the  low  Reynolds 
number  measurements,  but  at  high  angles,  the  theory  is  in  better 
agreement  with  the  high  Reynolds  number  results.  This  phenomenon 
is  not  completely  understood  at  this  time,  and  further  evaluation 
with  other  experimental  measurements  is  recommended. 
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The  effects  of  compressibility  on  the  measured  and  predicted 
normal  force  distributions  are  shown  in  Figure  16  for  a  =  15 
deqrees.  Predicted  compressibility  effects  exhibit  the  correct 
trend,  but  they  are  less  than  those  measured  on  the  nose.  The 
overshoot  at  the  shoulder  of  the  model  is  not  predicted. 


Measured  and  predicted  pressure  distributions  are  presented 
in  Fiqure  17.  The  increment  in  pressure  coefficient  due  to  com¬ 
pressibility  effects  appears  to  be  in  reasonable  agreement  with 
experiment,  although  there  is  some  disagreement  in  the  lower 
pressure  reqions  near  8  >  60  degrees  on  the  nose.  This  obviously 
accounts  for  the  problem  with  normal  force  discussed  above.  The 
decreasing  effect  of  compressibility  along  the  cylindrical  por¬ 
tion  of  the  body  is  in  good  agreement  with  experiment. 


The  above  comparisons  indicate  the  compressibility  correc¬ 
tion  included  in  the  prediction  method  and  code  represents  the 
measured  effects;  however,  further  verification  is  needed.  The 
limited  experimental  data  available  are  not  sufficient  to  define 
the  Mach  number  and  angle  of  attack  limits  of  the  code. 


Additional  comparisons  of  measured  and  predicted  pressure 
distributions  on  other  axisymmetric  bodies  are  presented  for 
further  verification.  For  these  comparisons,  the  2.0-  and  3.5- 
caliber  oqive-cylinder  models  tested  by  Lamont  (Refs.  17  and  24) 
over  a  wide  range  of  angles  of  attack  and  Reynolds  numbers  are 
used.  Only  experimental  data  at  moderate  anqles  of  attack  are 
considered  to  quarantee  flow  symmetry. 

In  Figure  18,  measured  and  predicted  circumferential  pres¬ 
sure  distributions  on  the  3.5-caliber  ogive-cylinder  model  at 
a  =  20  degrees  are  compared.  The  predictions  are  shown  with 

(solid  curve)  and  without  (dashed  curve)  vortex  induced  effects 
included.  As  in  the  previous  results,  the  separation-vortex  ef¬ 
fects  improve  the  shape  of  the  pressure  distribution  curves, 
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particularly  on  the  lee  side;  however,  there  is  a  general  lack  of 
agreement  between  the  experiment  and  theory.  There  appears  to  be 
a  uniform  shift  between  the  measured  and  predicted  pressure  coef¬ 
ficients,  even  on  the  windward  side  of  the  model  where  vortex 
induced  effects  should  be  negligible. 

In  Figure  19,  similar  comparisons  are  shown  for  the  2.0- 
caliber  ogive -cylinder  model  at  a  =  30  degrees.  Because  of  the 
higher  angle  of  attack,  vortex  induced  effects  are  larger  than 
before,  but  the  character  of  the  comparisons  is  nearly  the  same; 
that  is,  on  the  cylindrical  portion  of  the  model,  the  experiment 
and  theory  disagree  by  a  constant  increment  in  pressure.  There 
are  other  areas  of  disagreement  as  it  appears  that  the  predicted 
vortex  effects  are  often  stronger  than  exhibited  by  the  experi¬ 
ment. 


In  an  effort  to  better  understand  the  use  of  the  vortex 
cloud  prediction  method  and  in  an  attempt  to  resolve  the  differ¬ 
ences  discussed  above,  the  following  is  presented  for  discussion 
purposes.  Since  the  windward  meridian  should  have  minimal  vortex 
induced  effects  at  moderate  angles  of  attack,  a  comparison  of 
experiment  and  theory  along  this  meridian  should  be  a  measure  of 
the  accuracy  of  the  linear  or  attached-flow  part  of  the  theory. 
In  Figure  20,  a  number  of  experimental  results  for  various  ogive- 
cylinder  models  at  a  =  20  degrees  are  compared  with  a  single 
predicted  result.  The  ogive  noses  vary  from  1.5  to  3.5  diameters 
in  length  (Refs.  17,  20,  24,  and  25),  but  the  pressure  on  the 
cylinder  aft  of  the  nose  should  be  nearly  independent  of  the  nose 
shape.  The  predicted  pressures  are  for  the  3.5-caliber  ogive  of 
Lamont  and  the  3-caliber  ogive  of  Tinling  and  Allen.  It  is  ob¬ 
served  that  the  theoretical  results  are  in  agreement  with  the 
experiments  of  References  20  and  25  and  in  disagreement  with  the 
experiments  of  References  17  and  24.  The  authors  have  no  expla¬ 
nation  for  these  discrepancies. 
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At  higher  angles  of  attack,  the  separation  vorticity  becomes 
asymmetric  and  causes  both  an  induced  normal  force  and  side 
force.  As  discussed  in  a  previous  section,  an  asymmetric  solu¬ 
tion  develops  only  after  the  solution  is  perturbed,  and  the  final 
result  is  somewhat  dependent  on  the  perturbation.  Some  prelimi¬ 
nary  results  are  presented  to  illustrate  an  asymmetric  result. 

Detailed  flow  field  measurements  and  pressure  distributions 
on  a  3.0-caliber  ogive-cylinder  model  at  a  =  45  degrees  are 
available  in  Reference  26.  Based  on  earlier  experience  with 
asymmetric  solutions  (Ref.  2),  the  perturbation  region  was  lim¬ 
ited  to  the  first  diameter  of  length  of  the  nose,  and  the  magni¬ 
tude  of  the  perturbation  was  varied  to  find  an  appropriate  solu¬ 
tion.  In  Reference  2,  the  position  of  the  shed  vorticity  was 
used  as  the  criterion  to  evaluate  the  perturbation;  however,  in 
the  present  solution,  a  more  appropriate  criterion  is  the  fre- 
guency  of  variation  of  the  predicted  side-force  distribution.  In 
Figure  21,  the  measured  and  predicted  normal-force  and  side-force 
distributions  are  shown.  Two  different  perturbations,  A8  =  1  and 
2  degrees,  are  shown,  and  even  though  the  magnitude  of  the  pre¬ 
dicted  side  force  is  less  than  that  measured,  the  freguency  of 
the  side  force  matches  best  when  A8  =  2  degrees.  This  is  the 
perturbation  used  for  all  the  predicted  results  to  follow.  The 
predicted  normal  force  is  below  the  level  of  the  measured  normal 
force  by  a  small  amount,  but  the  side  force  is  considerably  less 
than  that  measured.  It  is  apparent  that  the  vortex  induced  ef¬ 
fects  are  too  weak  for  this  case. 

Measured  and  predicted  circumferential  pressure  distribu¬ 
tions  at  three  axial  stations  are  shown  in  Figure  22.  Since  the 
flow  is  asymmetric,  the  full  circumference  is  shown  in  these 
comparisons.  The  major  areas  of  disagreement  are  generally  for 
80°  <  8  <  160°  and  200°  <  8  <  270°;  although  at  the  aft  station 
where  the  error  in  side  force  is  greatest,  the  disagreement  is 
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concentrated  in  the  region  80°  <  0  <  160°  which  causes  the  maxi¬ 
mum  effect  on  predicted  side  force.  The  explanation  of  this 
problem  is  best  accomplished  by  a  discussion  of  the  separation 
wake  using  the  following  figures. 

The  predicted  vortex  cloud  at  x/D  =  4.7  is  shown,  in  Figure 
23  where  an  obvious  asymmetry  has  developed.  At  this  station, 
the  cloud  on  the  left  side  of  the  body  has  nearly  broken  away 
from  its  feeding  sheet  and  formed  a  free  vortex.  Reverse  flow 
vorticity  is  not  included  in  the  predicted  results  because  of 
numerical  difficulties  with  this  particular  case.  The  pressure 
distribution  on  the  lee  side  and  the  reverse  flow  velocity  field 
induced  by  the  primary  vortex  field  permits  the  calculation  of 
secondary  separation  and  the  shedding  of  secondary  vorticity. 
When  these  phenomena  were  included  in  the  calculation,  the  ef¬ 
fects  of  the  additional  vorticity  were  lost  when  the  primary 
vorticity  pulled  the  secondary  vortices  into  the  main  cloud  and 
effectively  eliminated  their  local  effect  on  the  pressure  distri¬ 
bution.  The  need  for  the  secondary  vorticity  is  illustrated  in 
the  following  figures. 

The  measured  velocity  field  at  the  same  axial  station  is 
shown  in  Figure  24(a)  where  the  individual  flow  vectors  are 
shown.  The  predicted  velocity  vectors  in  the  same  general  region 
of  the  lee  side  flow  field  are  shown  in  Figure  24(b).  Both  mea¬ 
sured  and  predicted  flow  fields  are  similar  in  magnitude  and  flow 
directions  except  in  the  reverse  flow  region  near  the  surface  of 
the  body.  It  is  obvious  in  Figure  24(a)  that  a  secondary  vortex 
field  exists  in  the  region  of  interest,  120°  <  0  <  170°.  The 
difficulties  with  the  predicted  secondary  separation  vorticity 
were  not  resolved  in  the  current  investigation. 

The  success  of  the  previous  flow  field  comparisons  stimu¬ 
lated  another  set  of  comparisons  to  be  sure  the  previous  results 
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were  not  fortuitous  nor  accidental.  Experimental  flow  field 
velocity  components  on  the  lee  side  of  a  3.5-caliber  ogive- 

cylinder  model  at  a  =  22.4  and  37.5  degrees  are  available  in 

Reference  27.  At  the  lower  angle  of  attack,  a  =  22.4  degrees, 

the  symmetric  flow  field  is  shown  in  Figure  25  for  three  dif¬ 

ferent  positions  in  and  near  the  vortex  cloud.  The  location  of 
the  experimental  results  and  the  predicted  vortex  cloud  are  shown 
in  the  sketches  in  the  figure.  In  each  case  the  predicted  down- 
wash  velocity  distribution  is  in  very  good  agreement  with  the 
experimental  results.  The  observed  roughness  when  the  predic¬ 
tions  are  carried  out  through  the  center  of  the  cloud  is  caused 
by  the  influence  of  a  discrete  vortex  too  near  the  calculated 
point. 

At  a  =  37.5  degrees,  the  flow  asymmetry  was  developed  with 
a  A0  =  2  degrees  perturbation  in  accordance  with  the  previous 
asymmetric  results.  The  velocity  field  comparisons  are  shown  in 
Figure  26,  and  the  location  of  the  results  and  the  asymmetric 
vortex  field  is  illustrated  in  the  sketches.  As  above,  the 
agreement  is  very  good.  The  roughness  of  the  predictions  is 
again  caused  by  the  near  proximity  of  several  discrete  vortices 
as  seen  in  the  sketches. 

Noncircular  Bodies 

The  noncircular  body  options  in  VTXCLD  have  been  tested  for 
a  number  of  different  cros:  section  shapes,  both  elliptic  and 
arbitrary,  to  verify  the  transformation  procedures.  Unfortu¬ 
nately,  only  a  limited  guantity  of  experimental  data  are  avail¬ 
able  for  comparison  purposes.  A  series  of  elliptic  bodies  and 
one  square  body  with  rounded  corners  are  examined  to  illustrate 
the  capabilities  of  the  code. 


Measured  normal-force  coefficients  on  a  series  of  elliptic 
cross  section  bodies  at  high  angles  of  attack  are  available  in 
Reference  21.  The  cross  section  shapes  are  1.5,  2,  and  3.5:1 
ellipses  with  egual  areas,  and  each  model  has  a  3-caliber  ogive 
nose  and  7-caliber  cylindrical  afterbody.  Measured  and  predicted 
normal-force  coefficients  are  shown  in  Figure  27.  Agreement 
between  experiment  and  theory  is  very  good  for  a  <  20  degrees  for 
all  models,  and  the  error  at  larger  angles  of  attack  appears  to 
be  nearly  the  same  for  the  three  cross  sections.  As  expected, 
the  nonlinear  effects  become  more  evident  as  the  ratio  of  major 
to  minor  axes  increases. 


The  predicted  vortex  cloud  patterns  at  the  end  of  the  body 
for  each  of  the  three  cross  sections  are  shown  in  Figure  28.  The 
dimensionless  strength  of  the  vorticity  on  the  positive  side  of 
each  body  is  also  shown  in  the  figure. 

Though  pressure  measurements  are  not  available,  it  is  inter¬ 
esting  to  observe  the  change  in  pressure  distribution  along  the 
elliptic  bodies  as  the  vortex  strength  increases.  In  Figure  29, 
the  predicted  pressure  distributions  around  the  2:1  elliptic  body 
at  a  =  20  degrees  for  three  axial  stations  are  shown.  The  solid 
curves  from  the  rational  flow  model  are  compared  with  the  poten¬ 
tial  flow  distributions,  and  it  appears  that  the  major  nonlinear 
effects  are  an  increase  in  the  minimum  pressure  at  the  side  of 
the  body  and  lower  pressure  on  the  lee  side  of  the  body. 

The  prediction  method  was  next  applied  to  a  square  body  with 
rounded  corners  at  a  =  20  degrees  and  M  =  0.5  from  Reference 

oc 

22.  The  model  has  a  3-caliber  nose  followed  by  a  10-caliber 
cylindrical  afterbody.  The  measured  and  predicted  normal-force 
coefficients  for  the  body  at  <p  =  0  and  45  degrees  roll  angles  are 
compared  in  Figure  30(a)  and  (b),  respectively.  At  <p  =  0°,  the 
predicted  normal  force  is  greater  than  the  measured  value.  It  is 
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interesting  that  the  measured  results  show  very  little  effect  of 
corner  radius;  that  is,  the  smallest  radius  tested  behaves  like  a 
circular  cross  section  at  subsonic  Mach  numbers. 

Measured  and  predicted  results  on  the  same  body  rolled  45 
degrees  are  shown  in  Figure  30(b).  The  agreement  is  very  good 
for  the  three  Mach  numbers  0.5,  0.7,  and  0.8.  In  the  rolled 

condition,  the  effect  of  corner  radius  is  significant  on  the 
measured  normal  force.  This  phenomena  may  be  explained  by  look¬ 
ing  at  the  vortex  shedding  results. 

In  Figure  31,  the  predicted  vortex  wake  on  the  square  body 
with  rounded  corners  is  shown  for  <j>  =  0  and  45  degrees.  At 

4>  =  o,  separation  can  occur  at  either  the  upper  or  lower  corner, 
and  it  may  occur  at  both  locations  simultaneously  if  the  lower 
portion  of  the  flow  reattaches  on  the  flat  side.  The  prediction 
method  cannot  handle  multiple  separation  or  reattachment;  there¬ 
fore,  it  is  not  surprising  that  the  normal-force  coefficients  in 
Figure  30(a)  are  in  poor  agreement.  When  the  model  is  rolled  45 
degrees,  there  is  only  one  adverse  pressure  region  near  the  cor¬ 
ner  on  the  side,  and  a  single  separation  location  is  indicated. 
The  prediction  method  appears  to  have  little  problem  with  this 
situation.  The  vortex  strengths  shown  on  each  figure  indicate 
that  the  vortex  wake  from  the  4>  =  45  degrees  flow  condition  is 
slightly  stronger  than  that  corresponding  to  the  unrolled  condi- 
t  i  c  ' . 


CONCLUSIONS 

An  engineering  prediction  method  based  on  a  rational  flow 
modeling  technique  and  the  associated  code  VTXCLD  to  predict  the 
vortex  shedding  from  circular  and  noncircular  bodies  in  subsonic 
flow  at  angles  of  attack  and  roll  are  described.  Comparisons  of 
measured  and  predicted  aerodynamic  characteristics  and  flow  field 
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quantities  are  used  to  verify  the  flow  model  and  prediction 
method  for  a  variety  of  conf iqurations  under  a  wide  range  of  flow 
conditions.  The  method  has  proved  successful  in  representing  the 
principal  features  of  the  complex  flow  field  on  the  lee  side  of 
missiles  at  moderate  angles  of  attack;  therefore,  it  has  applica¬ 
tion  as  an  engineering  or  preliminary  design  technique  directed 
at  the  prediction  of  nonlinear  aerodynamic  characteristics  re¬ 
sulting  from  high  angle  of  attack  flows. 

The  prediction  method  described  herein  has  further  applica¬ 
tion  as  one  component  of  a  larger  prediction  method  for  complete 
configurations  consisting  of  an  arbitrary  body  and  multiple 
fins.  The  ability  to  model  the  correct  flow  field  in  the  vicin¬ 
ity  of  the  body  leads  to  the  capability  to  calculate  body  separa¬ 
tion  wake  induced  interference  effects  on  fins  and  other  control 
surfaces.  The  vortex  shedding  analysis  and  prediction  techniques 
are  also  applicable  for  use  in  other  codes;  for  example,  the 
methods  developed  in  this  investigation  can  be  transferred  to 
higher  order  codes  such  as  panel  codes.  This  was  demonstrated  in 
the  supersonic  vortex  shedding  work  reported  in  Reference  3. 

Finally,  the  prediction  method  presented  in  this  report  must 
be  considered  as  preliminary.  Even  though  verification  of  the 
method  by  comparisons  with  experiment  has  demonstrated  some  suc¬ 
cess,  there  are  configurations  and  flow  conditions  which  must  be 
approached  with  caution.  For  this  reason,  several  recommenda¬ 
tions  for  improvements  to  the  method  are  suggested  in  the  follow¬ 
ing  section. 


RECOMMENDATIONS 


In  the  course  of  development  and  verification  of  the  code 
VTXCLD,  several  specific  areas  requiring  additional  work  were 
identified.  Since  the  additional  effort  was  beyond  the  scope  of 
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the  present  investigation,  they  are  noted  below  as  recommenda¬ 
tions  for  future  work. 

The  first  and  most  important  recommendation  is  for  a  thor¬ 
ough  testing  and  further  verification  of  VTXCLD.  The  code  should 
be  applied  to  a  variety  of  configurations  and  flow  conditions  for 
which  experimental  results  are  available  for  comparison  pur¬ 
poses.  This  should  better  define  the  operational  limits  of  the 
method.  The  major  difficulty  will  be  the  lack  of  useful  experi¬ 
mental  results  on  noncircular  bodies. 

As  noted  in  the  text,  use  of  the  code  at  high  angles  of 
incidence  for  asymmetric  vortex  wake  calculations  presents  a 
problem  regarding  the  required  perturbation  necessary  to  develop 
asymmetric  solutions.  It  is  recommended  that  the  code  be  applied 
to  a  number  of  asymmetric  cases  for  which  experimental  data  are 
available  in  an  effort  to  correlate  the  magnitude  and  location  of 
the  perturbation  in  light  of  comparisons  with  data.  This  is  the 
only  method  which  will  permit  the  code  to  be  used  with  confidence 
for  asymmetric  flow  conditions. 

A  similar  problem  exists  with  the  transition  option  included 
in  the  code.  Some  effort  should  be  expended  to  better  understand 
the  effects  of  changing  the  character  of  separation  from  laminar 
to  turbulent.  As  before,  correlation  of  predicted  results  with 
experiment  may  provide  the  necessary  rules  for  using  transition 
with  confidence. 

The  numerical  difficulty  associated  with  secondary  separa¬ 
tion  discussed  in  a  previous  section  needs  further  investiga¬ 
tion.  A  minor  modification  may  correct  the  problem,  but  some 
additional  effort  is  required  to  understand  the  source  of  the 
problem. 
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The  prediction  method  may  have  additional  uses  which  have 
not  been  explored.  For  example,  the  code  has  vortex  shedding 
technology  built  into  it  which  can  be  transferred  to  three-dimen¬ 
sional  panel  codes.  The  code  also  has  potential  use  in  providing 
separation  line  and  vortex  wake  information  for  starting  solu¬ 
tions  for  higher  order  prediction  methods. 

Nielsen  Engineering  &  Research,  Inc. 

Mountain  View,  CA  94043-2287 
November  1985 
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(b)  Numerical  transformation  procedure 
Figure  2.-  Conformal  mapping  nomenclature. 
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Sketch  of  crossflow  plane  separation  points. 


Figure  4.  -  Body  crossflow  plane  nomenclature. 
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Subroutine  calling  sequence  for 
program  VTXCLD. 
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Figure  6.-  Concluded 
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Figure  7.-  Common  block  cross  reference  list  for 
program  VTXCLD. 
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Figure  7 


Concluded 


Figure  8.-  Continued 
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Figure  8.-  Continued 


Plane  of  syranetry 


YRC 


Figure  gj.-  Convention  for  ordering  coordinates  for  a 
noncircular  cross  section  at  X  =  XFC (J) . 
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(a)  Sample  Case  1 

Figure  10.-  Sample  cases  for  program  VTXCLD. 
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Item 

(1)  0  0  1  1  0  1  1  0  7  0  0  0  0  0 

(2)  4  0  1  0  3  1  0  0  1  0 

(3)  SAMPLE  CASE  2 


NIELSEN  ENGINEERING  4  RESEARCH,  INC.  PROGRAM  VTXCLD 
3  CALIBER  OGIVE-CYLINDER  REF.  NASA  TN  0-1297 

INITIAL  FUM  CONDITIONS  ALPHAC  =  20.00  DEGREES,  PHI  *=  0  DEGREES 


(4) 

17.34900 

4.70000 

14.10000 

50.47800 

4.70000 

(5) 

20.00000 

0.00000 

440000. 

(7) 

2.35000 

47.000 

2.35000 

0.0 

0.0  1.05  0.0 

(8) 

.04 

0.0 

0.0 

0.0 

0.0  0.0  0.0 

(9) 

61  2 

(10) 

0.20000E-02  0.25126E-02  0.31565E-02  0.39648E-02  0.49785E-02  0.62473E-02 

0.78333E-02  0.98141E-02  0.12284E-01  0.15366E-01  0.19194E-01  0.23941E-01 
0.29802E-01  0.37014E-01  0.45830E-01  0.56547E-01  0.69476E-01  0.84914E-01 
0.10314E+00  0.12439E+00  0.14874E+G0  0.17615E+00  0.20641E+00  0.23895E+00 
0.27310E+00  0.30800E+00  0.34291E+00  0.37782E400  0.41273E+00  0.44765E+00 
0.48256E+00  0.51747E+00  0.55238E+00  0.58730E+00  0.62221E+00  0.65712E+00 
0.69203E+00  0.72695E+00  0.76186E+00  0.79677E+00  0.83168E+00  0.86660E+00 
0.90151ET00  0.93642E+00  0.97133E400  0.10062E+01  0.10411E+01  0.10753E+01 
0.11077E+01  0.11375E+01  0.11641E+01  0.11871E+01  0.12066E+01  0.12226E+01 
0.12356E+01  0.12459E+01  0.12539E+01  0.12601E-W1  0.12649E+01  0.12685E+01 
0.12712E+01 

(11)  0.81719E-07  0.26387E-07  0.11344E-06  0.13636E-06  0.24S35E-0S  0.3S439E-OG 
0.59532E-06  0.91349E-06  0.14535E-05  0.22591E-05  0.34840E-05  0.53899E-05 
0.81163E-05  0.12185E-04  0.17927E-04  0.25719E-04  0.35925E-04  0.48499E-04 
0.61993E-04  0.74917E-04  0.84224E-04  0.82314E-04  0.74762E-04  0.27384E-04 
0 .13014E-05-0 .35579E-04  0.29103E-04-0.32427E-04  0.303S4E-04-0.31699E-04 
0 .30799E-04-0 .31411E-04  0.31004E-04-0.31254E-04  0.31134E-04-0.31135E-04 
0 . 31252E-04-0 . 31 005E-04  0.31408E-04-0.30802E-04  0 . 31691E-04-0 . 30371E-04 
0.32373E-04-0.28832E-04  0.31G65E-04  0.13720E-04-0.44716E-04-0.77072E-04 
-0.92338E-04-0 .92578E-04-0 .81301E-04-0 .63716E-04-0 .46689E-04-0 .31532E-04 
-0 .20487E-04-0 . 12619E-04-0 .76726E-05-0 . 43306E-05-0 . 27866E-05-0 . 10763E-05 
-0.14198E-05 

(12)  71 


0.00100 

0.00119 

0.00141 

0.00168 

0.00199 

0.00237 

0.00281 

0.00396 

0.00470 

0.00558 

0.00662 

0.00786 

0.00932 

0.01104 

0.01548 

0.01831 

0.02163 

0.02554 

0.03011 

0.03545 

0.04166 

0.05721 

0.06679 

0.07775 

0.09019 

0.10424 

0.11995 

0.13737 

0.17726 

0.19953 

0.22310 

0.24769 

0.27298 

0.29857 

0.52902 

0.96991 

1.01552 

1.04105 

1.06611 

1.09032 

1.11331 

1.13480 

1.17248 

1.18849 

1.20261 

1.21492 

1.22554 

1.23462 

1.24231 

1.25419 

1.25870 

1.26243 

1.26551 

1.26804 

1.27012 

1.27183 

1.27436 

1.27529 

1.27604 

1.27666 

1.27716 

1.27757 

1.27790 

0.00834 

0.80041 

0.00048 

0.00057 

0.00068 

0.00081 

0.00096 

0.00135 

0.00160 

0.00189 

0.00224 

0.00265 

0.00313 

0.00370 

0.00514 

0.00605 

0.00710 

0.00831 

0.00971 

0.01130 

0.01311 

0.01742 

0.01992 

0.02263 

0.02553 

0.02857 

0.03168 

0.03478 

0.04049 

0.04285 

0.04472 

0.04597 

0.04653 

0.04655 

0.04655 

(b)  Sample  Case  2 


Figure  10.-  Continued. 


.6 

0.0 


0.00334 

0.01308 

0.U4888 

0.15650 

0.75946 

1.15457 

1.24878 

1.27322 

0.00114 

0.00437 

0.01515 

0.03775 

0.04655 
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(15) 


(30) 


0.04655 

0.04642 

0.04558 

0.04401 

0.02911 

0.02568 

0.02238 

0.01931 

0.00819 

0.00678 

0.00560 

0.00461 

0.00169 

0.00137 

0.00112 

0.00091 

0.34149 

0.34123 

0.34092 

0.34056 

0.33743 

0.33642 

0.33523 

0.33381 

0.32183 

0.31803 

0.31357 

0.30837 

0.26693 

0.25465 

0.24074 

0.22510 

0.11935 

0.09305 

0.06543 

0.03676 

0.00000 

-0.01802 

-0.04771 

-0.07699 

-0.20441 

-0.22429 

-0.24204 

-0.25770 

-0.30890 

-0.31491 

-0.31991 

-0.32406 

-0.33606 

-0.33733 

-0.33836 

-0.33920 

2.35 

4.70 

9.40 

21.150 

0.04181 

0.03907 

0.03594 

0.03257 

0.01650 

0.01398 

0.01177 

0.00964 

0.00378 

0.00310 

0.00253 

0.00207 

0.00074 

0.00060 

0.00049 

0.34013 

0.33962 

0.33901 

0.33829 

0.33213 

0.33015 

0.32782 

0.32507 

0.30230 

0.29525 

0.28710 

0.27770 

0.20763 

0.18830 

0.16711 

0.14409 

0.00737 

0.00000 

0.00000 

0.00000 

-0.10545 

-0.13272 

-0.15847 

-0.18243 

-0.27134 

-0.28312 

-0.29318 

-0.30171 

-0.32748 

-0.33030 

-0.33261 

-0.33451 

-0.33988 

-0.34044 

-0.34090 

28.20 

35.25 

47.00 

(b)  Concluded 
Figure  10.-  Continued. 
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Item 


0 


(1) 

1  0 

0  1 

0  1 

1  0 

8  0 

0  0 

0  0 

(2) 

5  0 

1  0 

3  1 

0  0 

1  0 

(3) 

SAMPLE  CASE  3 

NIELSEN  ENGINEERING  £  RESEARCH,  INC. 

PROGRAM  VTXCLD 

3  CALIBER  NOSE 

REF. 

HIT  AL  138 

2  :  1  ELLIPTIC  CROSS  SECTION 

INITIAL  FLOH  CONDITIONS  ALPHAC  «  15.00  DEGREES, 

PHI  =  0  DEGREES 

(4) 

4.4301 

2.375 

7.125 

23.75 

2.375 

(5) 

15.00000 

0.80000 

1348000. 

(7) 

1.1875 

23.750 

1.1875 

0.0 

0.0 

1.05 

0.0 

(8) 

.04 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

(9) 

-1  2 

(12) 

65 

(13) 

0.00000 

0.00023 

0.00092 

0.00208 

0.00369 

0.00576 

0.00829 

0.01468 

0.01854 

0.02284 

0.02756 

0.03270 

0.03825 

0.04421 

0.05729 

0.06440 

0.07188 

0.07970 

0.08787 

0.09636 

0.10517 

0.12366 

0.13333 

0.14325 

0.15341 

0.16380 

0.17440 

0.18520 

0.20730 

0.21857 

0.22997 

0.24147 

0.25307 

0.26474 

0.27646 

0.30000 

0.47500 

0.65000 

0.82500 

1.00000 

1.02354 

1.04693 

1.09271 

1.11480 

1.13620 

1.15675 

1.17634 

1.19483 

1.21213 

1.24271 

1.30000 

1.23579 

1.26730 

1.27716 

1.28532 

1.29171 

1.29631 

(14) 

0.00000 

0.00008 

0.00034 

0.00076 

0.00134 

0.00208 

0.00298 

0.00521 

0.00653 

0.00798 

0.00954 

0.01121 

0.01297 

0.01481 

0.01870 

0.02072 

0.02277 

0.02485 

0.02693 

0.02901 

0.03107 

0.03508 

0.03701 

0.03888 

0.04067 

0.04237 

0.04397 

0.04547 

0.04811 

0.04923 

0.05022 

0.05107 

0.05177 

0.05232 

0.05271 

0.05303 

0.05303 

0.05303 

0.05303 

0.05303 

0.05271 

0.05177 

0.14811 

0.04547 

0.04237 

0.03888 

0.03508 

0.03107 

0.02693 

0.01870 

0.00000 

0.01481 

0.01121 

0.00798 

0.00521 

0.00298 

0.00134 

(15) 

0.36495 

0.36463 

0.36367 

0.36209 

0.35987 

0.35703 

0.35358 

0.34490 

0.33969 

0.33393 

0.32762 

0.32080 

0.31347 

0.30566 

0.28868 

0.27955 

0.27003 

0.26013 

0.24987 

0.23928 

0.22838 

0.20573 

0.19402 

0.18207 

0.16992 

0.15756 

0.14503 

0.13234 

0.10654 

0.09346 

0.08029 

0.06703 

0.05371 

0.04033 

0.02691 

(16) 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

-0.02691 

-0.05371 

-0.10654 

-0.13234 

-0.15756 

-0.18207 

-0.20573 

-0.22838 

-0.24987 

-0.28868 

-0.36495 

-0.30566 

-0.32080 

-0.33393 

-0.34490 

-0.35358 

-0.35987 

0.00000 

0.00011 

0.00045 

0.00101 

0.00178 

0.00277 

0.00397 

0.00695 

0.00871 

0.01064 

0.01272 

0.01494 

0.01729 

0.01975 

1.02493 

0.02763 

0.03037 

0.03313 

0.03591 

0.03868 

0.04142 

0.04677 

0.04935 

0.05184 

0.05422 

0.05649 

0.05863 

0.06062 

0.06414 

0.06565 

0.06697 

0.06810 

1.06903 

0.06976 

0.07029 

0.07171 

0.07071 

0.07071 

0.07071 

8.07071 

0.07029 

0.06903 

0.06414 

0.06062 

0.05649 

0.05184 

0.04677 

0.04142 

0.03591 

0.02493 

0.01975 

0.01494 

0.01064 

0.00695 

0.00397 

0.00178 

( 

c)  Sample  Case  3 

Figure  10. 

Continued 

• 

.6 

0.0 


0.01126 

0.05056 

0.11427 

0.19617 

0.28822 

1.07003 

1.22812 

1.29908 

0.00402 

0.01673 

0.03309 

0.04685 

0.05295 

0.05022 

0.02277 

0.00034 

0.34954 

0.29739 

0.21719 

0.11950 

0.01346 

-0.08029 

-0.27003 

-0.36367 

0.00536 

0.02230 

0.04413 

0.06247 

0.07060 

0.06697 

0.03037 

0.00045 
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(17) 


(30) 


6.10000 

0.00000 

0.00006 

0.00022 

0.00050 

0.00347 

0.00435 

0.00532 

0.00636 

0.01247 

0.01381 

0.01518 

0.01657 

0.02339 

0.02467 

0.02592 

0.02711 

0.03207 

0.03282 

0.03348 

0.03405 

0.03535 

0.03535 

0.03535 

0.03535 

0.03207 

0.03031 

0.02825 

0.02592 

0.01247 

0.00987 

0.00747 

0.00532 

0.90000 

1.18750 

2.37SOO 

4.75000 

7.12500 

0.00089 

0.00139 

0.00199 

0.00268 

0.00747 

0.00864 

0.00987 

0.01115 

0.01795 

0.01934 

0.02071 

0.02206 

0.02825 

0.02931 

0.03031 

0.03123 

0.03451 

0.03488 

0.03514 

0.03530 

0.03535 

0.03514 

0.03451 

0.03348 

0.02339 

0.02071 

0.01795 

0.01518 

0.00347 

0.00199 

0.00089 

0.00022 

9.50000 

14.25000 

19.00000 

23.75000 

(c)  Concluded 
Figure  10.-  Continued. 
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Item 

(1)  2  9  1  1  0  1  1  0  8  0  0  0  0  0  1 

(2)  5110310811 

(3)  SAMPLE  CASE  4 

NIELSEN  ENGINEERING  i  RESEARCH,  INC.  PROGRAM  VTXCLD 

SQUARE  WITH  ROWED  CORNERS  REF.  SCKCIDER 

NUMERICAL  HAPPING 


INITIAL  aOH  CONDITIONS  ALPHAC  =  20.00  DEGREES,  PHI  *  0  DEGREES 


(4) 

12.5664 

4.00000 

12.00000 

52.00000 

4.8000 

(5) 

20.00000 

0.80000 

700000. 

(6) 

0.50000 

(7) 

2.00000 

52.000 

2.50000 

0.0 

0.0 

1.05 

0.0 

.6 

(8) 

.05 

0.0 

0.0 

0.0 

0.0 

0.8 

0.0 

0.0 

(9) 

-1  2 

(12) 

65 

(13) 

0.00000 

0.00018 

0.00071 

0.00160 

0.00284 

0.00443 

0.00638 

0.00866 

0.01129 

0.01426 

0.01757 

0.02120 

0.02515 

0.02942 

0.03401 

0.03889 

0.04407 

0.04954 

0.05529 

0.06131 

0.06759 

0.07412 

0.08090 

0.08790 

0.09513 

0.10256 

0.11019 

0.11801 

0.12600 

0.13416 

0.14246 

0.15090 

0.15946 

0.16813 

0.17690 

0.18575 

0.19467 

0.20365 

0.21266 

0.22171 

0.23077 

0.42308 

0.61538 

0.80769 

1.00000 

1.01811 

1.03610 

1.05387 

1.07131 

1.08831 

1.10477 

1.12058 

1.13564 

1.14987 

1.16318 

1.17548 

1.18670 

1.19676 

1.20562 

1.21320 

1.21947 

1.22439 

1.22793 

1.23006 

1.23077 

(14) 

0.00000 

0.00006 

0.00025 

0.00055 

0.00096 

0.00152 

0.00218 

0.00295 

0.00383 

0.00481 

0.00589 

0.00706 

8.00831 

0.00965 

0.01106 

0.01253 

0.01406 

0.01564 

0.01725 

0.01889 

0.02056 

0.02223 

0.02390 

0.02556 

0.02721 

0.02882 

0.03039 

0.03191 

0.03337 

0.03475 

0.03606 

0.03727 

0.03839 

0.03940 

0.04030 

0.04107 

0.04171 

0.04222 

0.04259 

0.04281 

0.04288 

0.04288 

0.04288 

0.04288 

0.04288 

0.04259 

0.04171 

0.04030 

0.03839 

0.03606 

0.03337 

0.03039 

0.02721 

0.02390 

0.02056 

0.01725 

0.01406 

0.01106 

0.00831 

0.00589 

0.00383 

9.00218 

0.00098 

0.00025 

0.00000 

(15) 

0.34584 

0.34522 

0.34460 

0.34357 

0.34213 

0.34027 

0.33800 

0.33530 

0.33218 

0.32864 

0.32467 

0.32026 

0.31542 

0.31014 

0.30442 

0.29825 

0.29162 

0.28455 

0.27701 

0.26902 

0.26056 

0.25165 

0.24227 

0.23243 

0.22212 

0.21136 

0.20014 

0.18847 

0.17635 

0.16379 

0.15081 

0.13740 

0.12357 

0.10936 

0.09475 

0.07978 

0.06445 

8.04880 

0.03283 

0.01657 

0.00000 

0.00000 

0.00000 

0.00000 

0.80000 

-0.03283 

-0.06445 

-0.09475 

-0.12357 

-0.15081 

-0.17635 

-0.20014 

-0.22212 

-0.24227 

-0.26056 

-0.27701 

-0.29162 

-0.30442 

-0.31542 

-0.32467 

-0.33218 

-0.33800 

-0.34213 

-0.34460 

-0.34584 

(19)  20  1 

(20)  26.00000 

(21)  25 

(22)  0.0800QE400  0.33333E400  8.66667E+00  O.lOOQOE+81  0.13333E+01  0.15885E+01 
0.18047E+01  0.19493E+01  0.20000E+01  0.20000E401  0.20000E+01  0.20000E+01 
0.20000E+01  0.28OOOE+O1  0.20000E+01  0.20000E+01  0.20000E401  0.19494E+01 
0.18047E401  O.15085E+O1  0.13333E+01  0.10000E401  0.66667E400  0.33333E+00 


(d)  Sample  Case  4 
Figure  10.-  Continued. 
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0.0W00E400 

( 23 )  -O.2OO00E4O1-O.2OOOOE+O1-O.200OOE+01-4.2O0OOE+OHI.2OO0OE+01-O.13493E+O1 

-0 .18047E401-0  .15885E401-0  .13333E+01-0  .1OOO0E+O1-O  .OO7E+OO-0 .  33333E+0Q 
O.00OOOE4M  O.33333E+0O  O.««7E+O0  0.10000E401  0.13333E-W1  0.13885E+01 
B.18047E+01  0.19493E401  0.20000E+01  8.20000E+01  Q.20000E+01  0.20000E+01 
0.20000E+01 

(30  )  2.00000  7.00000  12.00000  17.00000  22.00000  32.00000  42.00000  52.00000 


(d)  Concluded 
Figure  10.-  Continued- 
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Item 


(1) 

1  0 

0  1 

0  1 

1  0 

7  O 

0  0 

0  0 

0 

(2) 

4  • 

1  0 

3  1 

0  0 

1  1 

(3) 

SAMPLE  CASE  S 

COMPRESSIBILITY  CORRECTION 

NIELSEN  ENGINEERING  4  RESEARCH,  INC. 

PROGRAM  VTXCLD 

3  CALIBER  06KE-CYLINDER 

REF. 

NASA  TN  D-1297 

INITIAL 

FLOW  CONDITIONS  ALPHAC  *  15.00  DEGREES, 

PHI  *  0  DEGREES 

(4) 

17.34900 

4.70000 

14.10000 

50.47800 

4.70000 

(5) 

13.00000 

0.00000 

440000. 

(6) 

0.800 

(7) 

2.33000 

47.000 

2.35000 

0.0 

0.0 

1.05 

0.0 

.6 

(8) 

.04 

0.0 

0.0 

0.0 

8.0 

0.0 

0.0 

0.8 

(9) 

-1  2 

(12) 

71 

(13) 

0.00100 

0.00119 

0.00141 

0.00168 

0.00199 

0.00237 

0.00281 

0.00334 

0.00396 

0.00470 

0.00558 

0.00662 

0.00786 

0.00932 

0.01104 

0.01308 

0.01548 

0.01831 

0.02163 

0.02554 

0.03011 

0.03545 

0.04166 

0.04888 

0.03721 

0.06679 

0.07775 

0.09019 

0.10424 

0.11995 

0.13737 

0.15650 

0.17726 

0.19993 

0.22310 

0.24769 

0.27298 

0.29857 

0.52902 

0.75946 

0.98991 

1.01552 

1.04105 

1.06611 

1.09032 

1.11331 

1.13480 

1.15457 

1.17248 

1.18849 

1.20261 

1.21492 

1.22554 

1.23462 

1.24231 

1.24878 

1.23419 

1.25870 

1.26243 

1.26351 

1.26804 

1.27012 

1.27183 

1.27322 

1.27436 

1.27329 

1.27604 

1.27666 

1.27716 

1.27757 

1.27790 

(14) 

8.00034 

0.00041 

0.00048 

0.00057 

0.00068 

0.00081 

0.80096 

0.00114 

0.00133 

0.00160 

0.00189 

8.00224 

0.00265 

0.00313 

0.00370 

0.00437 

0.00314 

0.00603 

0.00710 

0.00831 

0.00971 

0.01130 

0.01311 

0.01515 

0.01742 

0.01992 

0.02263 

0.02353 

0.02857 

0.03168 

0.03478 

0.03775 

0.04049 

0.04285 

0.04472 

0.04597 

0.04653 

0.04655 

0.04655 

0.04655 

0.04655 

0.04642 

0.04358 

0.04401 

0.04181 

0.03907 

0.03594 

0.03257 

0.02911 

0.02368 

8.02238 

0.01931 

0.01650 

0.01396 

0.01177 

0.00964 

0.00819 

0.00678 

0.00560 

0.00461 

0.00378 

0.00310 

0.00253 

0.00207 

0.00169 

0.00137 

0.00112 

0.00091 

0.00074 

0.00060 

0.00049 

(15) 

0.34149 

0.34123 

8.34092 

0.34056 

0.34013 

0.33962 

0.33901 

0.33829 

0.33743 

0.33642 

0.33523 

0.33381 

0.33213 

0.33015 

0.32782 

0.32507 

0.32183 

0.31803 

0.31357 

0.30837 

0.30230 

0.29525 

0.28710 

0.27770 

0.26693 

0.25463 

0.24074 

0.22510 

0.20763 

0.18830 

0.16711 

0.14409 

0.11935 

0.09303 

0.06543 

0.03676 

0.00737 

0.00000 

0.00008 

0.80000 

0.00800 

-0.01802 

-1.84771 

-0.07699 

-0.10545 

-0.13272 

-0.15847 

-0.18243 

-0.28441 

-0.22429 

-0.24204 

-0.25770 

-0.27134 

-0.28312 

-0.29318 

-0.30171 

•0.38890 

-0.31491 

-0.31991 

-0.32406 

-0.32748 

-0.33030 

-0.33261 

-0.33451 

-0.33606 

-0.33733 

-0.33836 

-0.33920 

-0.33988 

-0.34044 

-0.34090 

(30) 

2.35 

4.70 

9.40 

21.150 

28.20 

33.23 

47.00 

(e)  Sample  Case  5 
Figure  10.-  Continued. 
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Item 


(1)  8  0  0  1  8  1  1  0  7  0  0  0  0  0  0 

(2)  4  8  1  0  3  1  0  0  1  0 

(3)  SAMPLE  CASE  6  TRANSITION  RESIGN 

NIELSEN  ENGINEERING  4  RESEARCH,  INC.  PROGRAM  VTXCLD 
3  CALIBER  OGIVE -CYLINDER  REF.  NASA  TN  D-1297 

INITIAL  FLOW  CONDITIONS  ALPHAC  -  20.00  DEGREES,  PHI  >  8  DEGREES 

(4)  17.34900  4.70080  14.10000  S0.47880  4.70000 

<5  )  20.00000  0.00800  440000. 

<7  )  2.35000  47.000  2.35000  7.05  11.75  1.05  0.0  .6 

<8)  .04  0.8  0.8  0.0  0.0  0.0  0.0  0.0 

(9)  61  2 

UO)  0.20000E-02  0.25126E-02  0.31565E-02  0.39648E-02  0.49785E-02  0.62473E-02 
0.78333E-02  0.98141E-02  0.12284E-01  0.15366E-01  0.19194E-01  0.23941E-01 
0.298Q2E-01  0.37014E-01  0.45830E-01  0.56547E-01  0.69476E-01  0.84914E-01 
0.10314E+08  0.12439E400  0.14874E+00  0.17615E+00  0.20641E+00  0.23895E400 
0.27310E+00  0.30800E+00  0.34291E+00  0.37782E+00  0.41273E+00  0.44765E+00 
0.48256E400  0.51747E+00  0.55238E+00  0.58730E+00  0.62221E+00  0.65712E+00 
0.69203E+08  0.72695E+00  0.76186E+00  Q.79677E+00  0.83168E+00  0.86660E+00 
0 .90151E+08  0.93642E+00  C.97133E+00  0.10Q62E+01  0.10411E+01  0.10753E+01 
0.11077E401  0.11375E+01  A.11641E+01  0.11871E+01  0.12066E+01  0.12226E+01 
0.12356E+01  0.12459E+01  0.12539E401  0.12601E401  0.12649E+01  0.12685E401 
0.12712E401 

(11)  0.81719E-87  0.26387E-07  0.11344E-06  0.13696E-06  0.24535E-06  Q.36439E-06 
0.59532E-06  0.91349E-06  0.1453SE-05  0.22591E-05  0.34840E-05  0.53899E-05 
0.81163E-05  0.12185E-04  0.17927E-04  0.25719E-04  0.35925E-04  0.48499E-04 
0.61993E-84  0.74917E-04  0.84224E-04  0.82314E-04  0.74762E-04  0.27384E-04 
0.13814E-05-0.35579E-04  0.29103E-04-0.32427E-04  0.30364E-04-0.31699E-04 
0.30799E-04-0.31411E-04  0.31004E-04-0.31254E-04  0.31134E-04-0.31135E-04 
0 . 31252E-04-0 .3100SE-04  0.31400E-04-0.30802E-04  0.31G91E-04-0.30371E-04 
0.32373E-04-0.28832E-04  0.31665E-04  0.13720E-04-0.4471GE-04-0.77072E-04 
-0 .92338E-04-0 . 92570E-04-0 . 81 301E-04-0 . 63716E-04-0 . 46689E-04-0 .31532E-04 
-0 . 20487E-04-0 . 12619E-04-0 . 76726E-05-0 . 43306E-05-0 . 27866E-O5-0 . 10763E-05 
-0.14198E-05 

(12)  71 

(13  )  0.00180  0.00119  0.00141  0.00188  0.00199  0.00237  0.00281  0.00334 

0.00396  0.00470  0.00558  0.00662  0.00786  0.00932  0.01104  0.01308 

0.01548  0.01831  0.02163  0.02554  0.03011  0.03545  0.04166  0.04888 

0.05721  0.06679  0.07775  0.09019  0.10424  0.11995  0.13737  0.15650 

0.17726  0.19953  0.22318  0.24769  0.27298  0.29857  0.52902  0.75946 

0.98991  1.01552  1.04105  1.06611  1.89032  1.11331  1.13480  1.15457 

1.17248  1.18849  1.20261  1.21492  1.22554  1.23462  1.24231  1.24878 

1.25419  1 .25670  1.26243  1.26551  1.26804  1.27012  1.27183  1.27322 

1.27436  1.27529  1.27604  1.27666  1.27716  1.27757  1.27790 

(14  )  0.80034  0.00041  8.00048  0.00057  0.00068  0.00081  0.00096  0.00114 

0.00135  0.00160  0.80189  0.00224  0.00265  0.00313  0.00370  0.00437 

0.00514  0.00605  0.00710  0.00831  0.00971  0.01130  0.01311  0.01515 

0.01742  0.01992  0.02263  0.02553  0.02857  0.03168  0.83478  0.03775 

0.04049  0.04285  0.04472  0.04597  0.04653  8.04655  0.04655  0.04655 

(f)  Sample  Case  6 
Figure  10.-  Continued. 
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(15) 


(30) 


0.04655 

0.04042 

0.04558 

0.04401 

1.02911 

0.02568 

0.02238 

0.01931 

0.10819 

0.00678 

0.00560 

0.00461 

0.00109 

0.00137 

0.00112 

0.00091 

0.34149 

0.34123 

0.34092 

0.34056 

0.33743 

0.33642 

0.33523 

0.33381 

0.32183 

0.31803 

0.31357 

0.30837 

0.20693 

0.25465 

0.24074 

0.22510 

0.11935 

0.09305 

0.06543 

0.03676 

0.00000 

-0.01802 

-0.04771 

-0.07699 

-0.20441 

-0.22429 

-0.24204 

-0.25770 

-0.30890 

-0.31491 

-0.31991 

-0.32406 

-0.33006 

-0.33733 

-0.33836 

-0.33920 

2.35 

4.70 

9.40 

21.150 

0.04181 

0.03907 

0.03594 

0.03257 

0.01650 

0.01398 

0.01177 

0.00984 

0.00378 

0.00310 

0.00253 

0.00207 

0.00074 

0.00060 

0.00049 

0.34013 

0.33962 

0.33901 

0.33829 

0.33213 

0.33015 

0.32782 

0.32507 

0.30230 

0.29525 

0.28710 

0.27770 

0.20763 

0.18830 

0.16711 

0.14409 

0.00737 

0.00000 

0.00000 

0.00000 

-0.10545 

-0.13272 

-0.15847 

-0.18243 

-0.27134 

-0.28312 

-0.29318 

-0.30171 

-0.32748 

-0.33030 

-0.33261 

-0.33451 

-0.33988 

-0.34044 

-0.34090 

28.20 

35.25 

47.00 

(f)  Concluded 


Figure  10.-  Continued. 
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Item 

(1) 

•  0 

1  1 

0  1 

1  0 

7  57 

0  1 

0  0 

1 

(2) 

5  1 

1  1 

3  1 

f  0 

0  0 

(3) 

WflE  CASE  7  ASttfCTMC  SEPARATION 

NIELSEN  ENGINEERING  &  RESEARCH,  INC. 

PROGRAM  VTXOJD 

3  CALIBER  06M-CYIIHCR 

REF. 

NSC  TR82-384 

ASWfCTRIC  SEPARATION  CASE 

INITIAL  R(W  COCITIQNS  ALPHAC  >  45.0  DEGREES,  1 

PHI  >  0  DEGREES 

(4) 

25.65207 

5.71500 

17.14500 

72.00900 

5.71500 

(5) 

45.60000 

0.00000 

106100. 

(  7) 

0.85725 

32.5735 

0.85725 

0.0 

0.0 

1.05 

0.0 

0.6 

( 8) 

.64 

0.0 

0.0 

5.145 

2.00 

0.0 

0.0 

0.0 

(  9) 

-1  2 

(12) 

55 

(13) 

0.00000 

0.00033 

0.00130 

0.00293 

0.00520 

0.00811 

0.01165 

0.01581 

0.02058 

0.02595 

0.03190 

0.03841 

0.04547 

0.05306 

0.06116 

0.06974 

0.67878 

0.08826 

0.09815 

0.10842 

0.11905 

0.13000 

0.14125 

0.15277 

0.16452 

0.17647 

0.18859 

0.20085 

0.21321 

0.22563 

0.23810 

0.48810 

0.73810 

0.98810 

1.23810 

1.25678 

1.27534 

1.29368 

1.31167 

1.32921 

1.34619 

1.36250 

1.37804 

1.39273 

1.40645 

1.41914 

1.43072 

1.44110 

1.45024 

1.45807 

1.46454 

1.46961 

1.47326 

1.47546 

1.47619 

(14) 

0.00000 

0.00011 

0.00045 

0.00100 

0.00176 

0.00273 

0.00389 

0.00522 

0.00672 

0.00836 

0.01013 

0.01200 

0.01396 

0.01598 

0.01804 

0.02012 

0.02219 

0.02423 

0.02622 

0.02814 

0.02997 

0.03168 

0.03327 

0.03471 

0.63599 

0.03709 

0.03801 

0.03874 

0.03926 

0.03958 

0.03968 

0.03968 

0.03968 

0.03968 

0.03968 

0.03944 

0.03874 

0.03758 

0.03599 

0.03401 

0.63168 

0.02907 

0.02622 

0.02321 

0.02012 

0.01701 

0.013% 

0.01106 

0.00836 

0.00595 

0.00389 

0.00222 

0.00100 

0.00025 

0.00000 

(15) 

0.34286 

0.34233 

0.34076 

0.33815 

0.33452 

0.32988 

0.32426 

0.31770 

0.31021 

0.30185 

0.29265 

0.28266 

0.27191 

0.26046 

0.24834 

0.23561 

0.22231 

0.20849 

0.19419 

0.17946 

0.16434 

0.14886 

0.13308 

0.11702 

0.10073 

0.08424 

0.06758 

0.05080 

0.03392 

0.01696 

0.00000 

0.00000 

0.00000 

0.60000 

0.00000 

-0.02545 

-0.05080 

-4.07593 

-0.10073 

-0.12508 

-4.14886 

-0.17195 

-4.19419 

-4.21547 

-4.23561 

-0.25448 

-0.27191 

-4.28775 

-0.30185 

-4.31407 

-4.32426 

-0.33232 

-0.33815 

-0.34168 

-0.34286 

(30) 

7.4295 

14.859 

20.574 

24.0030 

26.8605 

30.00375 

32.5755 

(31) 

2.12159 

2.12159 

1.92860 

2.29842 

1.72094 

2.45776 

1.50019 

2.59840 

1.76801 

2.71926 

1.02619 

2.81943 

0.77655 

2.89814 

0.52101 

2.95479 

0.26150 

2.98896 

0.00000 

3.00038 

-0.26150 

2.98896 

•4.32101 

2.95479 

-0.77655 

2.89814 

(g)  Sample  Case  7 
Figure  10.-  Continued. 
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-1.12619 

-1.26801 

-1.50019 

-1.72094 

-1.92860 

-2.1358 

2.62672 

2.38780 

2.13069 

1.85738 

1.56992 

1.27052 

0.96145 

0.64506 

0.32376 

0.00000 

-0.32376 

-0.64506 

-0.96145 

-1.27052 

-1.56992 

-1.85737 

-2.13069 

-2.38779 

-2.62672 

3.13186 

2.84699 

2.54044 

2.21456 

1.87183 

1.51485 

1.14634 

0.76911 

0.38602 

0.00000 

-0.38602 

-0.76911 

-1.14634 

-1.51485 

-1.87183 

-2.21456 

-2.54044 

-2.84699 

-3.13186 


2.81943 

2.71926 

2.59640 

2.45776 

2.29842 

2.12159 

2.62672 
2.84566 
3.04294 
3.21707 
3.36671 
3.49072 
3.58817 
3.65831 
3.70061 
3.71475 
3.70061 
3.65831 
3.58817 
3.49072 
3.36671 
3.21707 
3.04295 
2.84566 

2.62673 

3.13186 
3.39291 
3.62813 

3.83573 
4.01415 
4.16202 
4.27821 
4.36184 
4.41227 
4.42912 
4.41227 
4.36184 
4.27821 
4.16202 
4.01415 

3.83574 
3.62813 
3.39291 

3.13187 


(g)  Concluded 
Figure  10.-  Continued. 
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Item 


(1)  2Q0101100000000 

(2)  5  1  1  1  3  2  0  0  0  0 

(3)  SAtflE  CASE  8 

NIELSEN  ENGINEERING  £  RESEARCH,  INC.  PROGRAM  VTXCU) 

LOBE  BODY 
NUMERICAL  HAPPING 

INITIAL  FLQH  CONDITIONS  ALPHAC  =  15.00  DEGREES,  PHI  =  0  DEGREES 

(4)  15.05941  4.37884  7.00000  14.00000  4.37884 

(5)  15.00000  0.00000  729000. 

(7)  2.00000  14.000  2.00000  0.0  0.0  1.05  0.0  1.0 

(8)  .os  o.O  0.0  0.8  0.0  0.0  0.0  0.0 

(9)  o  0 

(12)  2 

(13)  o.ooooo  1.00000 

(14  )  0.15639  0.15639 

(I5)  0.00000  0.00000 

(19)  20  1 

(20  )  2.00000 

(21)  28 

(22)  O.OOOOOETOO  9.01558E+00  0.03105E+00  0.04629E+00  0.06110E+00  0.07524E+00 
0.08840E+00  0.10026E400  0.11049E400  0.11882E+00  0.12508E+00  0.12926E+00 
0.13147E+00  0.13194E+00  0.13097E+00  0.12891E+00  0.12607E+00  0.11400E+00 
0.10000E+00  0.08365E+00  0.07847E+00  0.07212E+00  0.06431E+00  0.05476E+00 
0.04333E+00  0.03011E+00  0.01546E+00  0.80000E+00 

(23)  HU7874E4OO^.178O9EMOHU7€12EM0HU7275E«<HU6786E«O-O.16134E+OO 
-0 .15312E+00-0 .14319E+00-0 .13168E+00-0 .11882E+00-0 .10496E+00-0 .09051E+00 
-0 .07590E400-0 . 06152E+OO-O .04767E400-0 .03454E+00-0 .02223E+00  0.0190QE+00 

0.05800E+00  0.09970E+00  0.11206ET00  0.12492E+00  0.13791E+00  0.15044E+00 
0.16171E+00  0.17077E400  0.17668E400  0.17874E+00 


(h)  Sample  Case  8 


Figure  10.-  Continued. 


-117- 


Item 


(1)  000181104090900 

(2)  5  0  1  0  3  1  0  0  1  0 

(3)  SAMPLE  CASE  9 

NIELSEN  ENGINEERING  4  RESEARCH,  INC.  PROGRAM  VTXCLD 

3  CALIBER  OGIVE-CYLINDER  REF.  NASA  TN  D-1297 

RESTART  CALCULATION 

INITIAL  FLON  CONDITIONS  ALPHAC  =  20.00  DEGREES,  PHI  =  0  DE6REES 


(4) 

17.34900 

4.70000 

14.10000 

50.47800  4.70000 

(5) 

20.00000 

0.00000 

440000. 

(7) 

21.15000 

47.000 

2.35000 

0.0  0.0 

1.05 

0.0  .6 

(8) 

.04 

0.0 

0.0 

0.0  0.0 

0.0 

0.0  0.0 

(9) 

61  2 

(10) 

0.20C30E-02  0.25126E-02  0.31565E-02  t  39648E-02  0.49785E-02  0.62473E-02 

0.78333E-02  0.98141E-02  0.12284E-01  0.15366E-01  0.19194E-01  0.23941E-01 
0.29802E-01  0.37014E-01  0.45830E-01  0.56547E-01  0.69476E-01  0.84914E-01 
Q.10314E+00  0.12439E+00  0.14874E+00  0.17615E+00  0.20641E+00  0.23895E+0Q 
0.27310E+00  0.30800E+00  0.34291E+00  0.37782E+00  0.41273E+00  0.44765E+00 
0.4825^+00  0.51747E+00  0.55238E+00  0.58730E+00  0.62221E+00  Q.65712E+00 
0.69203E+00  0.72695E+00  0.76186E+00  0.79677E+00  0.83168E+00  0.86660E+00 
0.90151E+00  0.93642E+0Q  0.97133E+00  0.10062E+01  0.104UE+01  0.10753E+01 
0.11077E+01  0.1137SE+01  0.U641E+01  0.11871E+01  0.12066E+01  0.12226E+01 
0.12356E+01  0.12459E+01  0.12539E+01  0.12601E+01  0.12649E+01  Q.12685E+01 
0.12712E+01 

( 11 )  0.81719E-07  0.26387E-07  0.11344E-O6  0.13696E-06  0.24535E-06  0  36439E-Q6 

0.59532E-06  0.91349E-06  0.14535E-05  0.22591E-05  0.34840E-05  0.53899E-05 
0.81163E-05  0.12185E-04  0.17927E-04  0.25719E-04  0.35925E-04  0.48499E-04 
0.61993E-04  0.74917E-04  II.84224E-04  0.82314E-04  0.74762E-04  0.27384E-04 
0 . 13014E-05-0 . 35579E-04  0.29103E-04-0.32427E-04  0.30364E-04-0 .31699E-04 
0.30799E-04-0.31411E-04  0.31004E-04-0.31254E-04  0.31134E-04-0.31135E-04 
0.31252E-04-0.31005E-04  0.31408E-04-0.30802E-04  0.31691E-04-0.30371E-04 
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Figure  10.-  Concluded. 
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INPUT  GEOMETRY  CHARACTERISTICS 
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Figure  11.-  Continued. 
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Figure  11.-  Continued 


80URCE  SPACING  18  0.79000  TIMES  LOCAL  RADIUS 
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Figure  11.-  Continued. 
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Figure  11.-  Continued. 


VELOCITIES  ON  THE  BOOT  SURFACE 
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Figure  11.-  Continued. 
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Figure  11.-  Continued. 
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Figure  11.-  Continued 


800V  SURFACE  PRESSURE  DISTRIBUTION 


■o 

>  RT 

q.  n 
o  . 
x  O 

I 


z  — 
a.  — 
u  . 
x  O 


r> 
O 
DC  ♦ 
H  Ui 
Ui  O 
«  >0 
r> 


»-NMf1«(Dir--OOift'rO,«'^ONMftO«NNO'0^O’M>OO^^OHinNN 

QOlM^nirnOflCaO^OO'OtttD'rO^lDH'OtQOHOtDiftOPl’tnH-OO-O 

MfmNNHftnniN-^w^TOOOOOOOOOH-^HHHriNWWWNfiM 


aooo  o'  odd  o'  boob  00666  606660  60000000000000 
a  r  1  1  i  1  t  >  1  1  1  1  i  1  1  1  1  1  1 


nrtoao^h«ro'^cD[DO^Ort«o«ta)o-ort(Dn«i>rtif)f»N'OrJO^rt 
(►'Orv-»-»«<oran-»o*cDr^(D-^far^ch'RO^-RrooiriiO‘r)r'Ch-»n*oO'Onr«i 
CL^r^'0|°nor>RrO'0r4rs'na'riry^'r  -o  o-  -*  r<  n  c  •rnnni'i-oO'iMDma) 

ORr^VRr^^pinr>wr4'-*-^ooooooo*R*^*R-R^-^^^-»-R»R-ROOOoo 


bbbbbbbbddbbdbbbdbbbbbbbboooooooooooo 

1  1  1  1  1  1  1  1  1  t  1  1  1  1  f  <  1  1  1  1  1 


OOOOftN-niO-ON«-‘t(Drt4«'On'ON^ttO'r-««NnO‘f4^N«000 

>-*-om-o<M)4W-<o^ajO‘NOft40>^«MrftOhO‘a>^a)^ft4>oO‘-oo-o-< 

s44nNff-rtfto-n-0(DON«’n«0<0N*0>0«<rtOD'0n-<(Drtn(0'M04-^ 

Knrtw«f4nnn<4ttn«ownif)«ni<inoftD4^^4nnnnw«ttr4 

^dddddddddddddo dodo dodo od do ddddddo dodo 


ocDO^o-hNOeO'dajn-<joo<to*o<HnO'0(D'r''^«ffi^'rffl-oonooa) 

>-»Nr4  0«ooorvcor^©©r^n*roor4fvRr'r«o*f^&*^o-o*o*a)©-*rvCDORr-t 
N4-nr4  0>4)r<(8M«oHro-o«NW«oO‘-in<^nn»oa)«ni^^n;;m««f5 
3rtrt«HHHOOOOHWrtnn44»rin«ftnft«  *v<  «Mfnnnnrtttmcvw 

dddddddddddddddddddddddddoddddddddddd 
1  1  1  1  1  1  1  1 

OOOfiO-OH(M1ffl^n'flt<OON(HDO-OHO‘^N^<l(D4l)ni^NO«0;0 

>o<JOM0‘Ofl*«O0'O-*nNnn'0fi-'00'4OO«r«&'N0*4n'Cn'0-0no 

\0'OnO‘404tD-«W4tn-«(h<ON04(h40'OriH45(DO'00^5jBtNO 

>ooHHNnnnf4t4«4nnnon--«HoooooooH^oooooo 

6666066606666  do66o6ddoddo<j>oocioo<ji. 


'^(D(DCD(D(DQ0(D(D(0(0a>(D(DQ0(0Q0CDC3a)09C0Ot0O000BCDQ3(DS>(DC003(DC000Ca 

dddddddddddddddddoddddddddddddddddodd 


0000000000000000000000000000000000000 

<0000000000000000000000000000000000000 

1-0000000000000000000000000000000000000 

. . 

O000*>0n000n0rton0fi0«0^0fl0«0fl0rt0ft0ft0ft0n0 

-HwtKnn4t»)«^'ONNOfflj-o*oo*<*«HNnn44fin-0'ONMD 


-flo«M>o4ntoft<BWNnoooflonftritD<nnj0  4o»N504 
HBo-ONnna-oNOO'flj'Oo-KnoON^ofrftfrONOj-nnMOfflH 
nw«oo<on»*>o*RNi«co»®ci-oo4«®ftO'0-NiJin'S®pwwn 
r.KNN'0'04>tf>tf)«)<r^nnrti'R^«ooo^^cinn<CRrr>nn<d'0<orvNr^rs. 


0000c 
1  1  1  T  i 


^?????????????OOQOOOOOOOOOOOOOOO° 


8oonMN04nn4n4OANAo<ooANAO4n4nn4ONNnpOi 
nh>O>«>ONOBOON-0OB*«BO-0Nnn^ONOBB>O0;h*n< 

O'0ciorf*O'0~r-‘-'0frn'0tDoc4r»n(¥C4©ffl-0n£5~*N 

OOHHfinnttBBBB'O'ONNNKhNh-O'OBfinBt 


h4)Q  ABB' 

<rnnw-»-( 


6660666666660006660666660000000000000 


•  r«nRrtf»4iiN>O0'O 


O-*rtO«rO-flND>O-fHn?«'0NOfrO:;2njJJ4N 

MHw-^«H*.HHrtwrt(KfiiNrtf4fiMKnnnnn(inn 


o 

o 

•f 

Ui 

o 
J  o 
©  o 
u  . 


O' 

O 

I 

UI 
© 
m 
tr  th 
©  . 


rx 

o 

1 

w 

O' 


rx 

o 

1 

Ui 

N 


O' 

O 

I 

u 
«-• 
OD 
>  r> 
o  . 
(4 


Z 

o 

►-* 

H 

< 

cr 

o 

ut 

»- 

z 


C<i 

O 

I 

UI 

* 

U 
Z  O' 
o  . 
W 


« 

o 

1 

-»  UI 
X  ^ 
SI 
<  4 
O  • 
ft 


u. 

u. 

UI 

o 

o 

H 

Z 

UI 


0* 

o 

I 

~  UI 
x  4 
^  ft 

>  o 
o  . 

*■ 


w 

? 

~  III 

x  n 
~  0^ 

X  ft 
o  . 

* 


ft 

s 


UI 

u 

<r 

o 

u. 


131 


STRATFORD  SEPARATION  CRITERION  (LAMINAR)  F<3)  -0.  022D2 


X  X 


Q 

o  o 

o 

o  o 

\ 

o  o 

s 

o  o 

■» 

o  o 

•» 

o  o 

a. 

CL 

o 

o  o 

O 

d  6 

a 

a 

■* 

o 

o 

o 

o 

CL 

c 

o 

0- 

o 

o 

u 

o 

o 

u 

o 

o 

o 

d 

d 

o 

tr> 

n 

C4 

o 

n 

o 

CM 

o 

o 

CL 

* 

o 

a 

CD 

o 

o 

O 

■«r 

o 

o 

O 

o 

o 

■ 

6  6 

o 

o 

* 

o 

o 

X 

o 

O  Cl  o 

0  0-0 

O 

a: 

O'  n  o 

C  O'  o 

o 

< 

-  <o  o 

<  —  o 

o 

> 

a 

0-^0 

o  o 

o 

o 

a 

z 

a 

< 

o  o  o 

<  o  o 

o 

cn 

►— 

o  o  o 

*-  o  o 

o 

h- 

in 

o  o  o 

UIOO 

o 

Z 

a 

a  .  . 

odd 

o  o 

o 

o 

01 

CD 

£ 

O 

n 

z 

n 

o 

CM  'O  o 

CM  O 

o 

n i 

h- 

n-oo 

n  o 

o 

< 

r*  n  o 

**  r*  o 

o 

cr 

• 

< 

o  o  o 

o  o 

d 

a. 

i 

X 

Oi 

cn 

< 

o 

z 

0*0 

o  o 

o 

a 

< 

o  n  o 

o  o 

o 

z 

> 

O  'O  o 

>  o  o 

o 

z 

a 

o 

odd 

o  o 

o 

u. 

H 

►- 

D 

o 

a 

z 

cr 

z 

h-  Ul 

1-  Ui 

o 

01 

cl  a: 

a.  <r 

>“• 

D 

D 

»“ 

o 

Z  U)  z 

Z  0) 

z 

< 

0  0)0 

o  m 

o 

ec 

Ul 

mU)m 

M  U| 

< 

cr 

h  a:  *•* 

$-  cr 

I- 

a 

D 

<  0.  < 

<  Q. 

< 

Ul 

cn 

Z  a: 

Z 

K 

CD 

(0 

o  .  < 

o  . 

< 

Ul 

<  Z  CL 

<  z 

CL 

a: 

H  ~  UJ 

UJ 

CL 

cn  E  CD 

cd  r 

(0 

U. 

o 

ui 

ui 

a 

a 

> 

a: 

m 

CD 

< 

r 

> 

> 

E 

* 

1 

o 

cn 

-JL32- 


Figure  11.-  Continued. 


BODY  SURFACE  PRESSURE  DISTRIBUTION 


«-krio^(>'-.NcO'<rr'‘-o<ror'--or^nr»non^nr'v-or'-0'*>ON’<rcDr^~.fr--»-or> 
Q^O(»ncDOHr.O'^NrsonnMDO-o(MDMnnONfj'OC-'<Hoa)na)0** 
\rar«i  —  -«ooo(D’Oio»rrrf-^0'N»r»n-.o-*ninN^-»rj^in>ooooo-*'-«r4rj 
Hfirtn»nirnc<--HH--.-ooooooooooo-^-4--H«rjriNNHN 

X . 

CLOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
Q  I  1  1  \  I  I  »  I  I  I  I  I  I  I  I  »  1  I 


-«rooo'f-CDCOO<a''0  —  fyr>-ioor^<r-oora»rn«tr^r-^-0'0<<r-ocD'Or>c>'rjh-io 
inra^‘<hO'Ocor»ci'00‘'*-»nc4'OCDNr>iRr-*«rncDa“r‘r<m^)-ooRr^«f>K«tDr^ 
c,-<'<ocoN'r^atf>*Th^on-oo*fjoM>0'’^'-'-*Hoa'iiN-or)vnnrtrii 
unnnf<f<w«*«^^oooooo-^^«fiirgrtNrdrliHTM^HH^Hw^- 

6666066606006606660 o 066066  odooddddbod 

t  I  I  I  I  I  I  I  I  !  I  S  I  I  I  t  I  !  I  I  I  t  I  I 

or^'0«ntnRrn>ooir<(>r.nr4n-.r^r^rjordr^r^— nrwnr^o-rani-onvmm-oh. 
>mo^n^oo^^(h-o-nni(DO(MONooomrjnH^O‘w^oo«invor) 
\o-rt«fNon'00‘««tM)'^H^vooio<<ni-0‘N<r^0‘4noN'r«^o 
Krt«Hf*mnnnn<r»<^««ioif)iftnniftin«n<^^^onnnNNNN« 

0606666666666660666666666666666666060 


ohOOOoo(D(Dnonn^«NNntD-o>fla)N^H^nonif>CD<j(Hh^Na‘N 

n  -  n  >o  ^coNconnffiOn 

\oom-ono‘^o«HhwcDnorH>0‘Hnnp)«HO.«oooNnoN«n-ioo 
3wm**^-«oooOH-.wwnn'>tt«ionnnrtv«r<firnnnNttNNw« 

ooooooo666dd6odddo666d666ddd6dddd6666 

ooo-n*r^p^rxon<-*ninrNnr^rs'«-*-^a-rNO-^nocDorvrjn>ooinnowcDOO 

>oo<^rni-o«nh'-*rvRrc^-^r<<ro'Co-*oooRrO'N«-*~o<rKnr<r>^^o<iwNO 

\o-or»CDna)Oi>o(DO^-oCD^rjo^O'0-JDN^^Mh-cnrtni*<(hNOMo 

>oo-^Nnnnn^'f'r«rnnnnfjN-Hoooooo-^-HHOoooo 

6666666666666666666666666666666666666 

1  1  1  1  1  1  1  1  1  1  1  1  1 

rT’*  +  'T-*'t*t*t*t*t* 

>  rtf  01  01  ot  01  c*  ex  01  o<  ex  ex  rn  oi  01  oi  01  01  ct  ex  ot  01  ex  ex  ex  ex  ni  oj  01  ni  oi  ni  ot  oi  01  m  ex  rtf 

D . 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


0  . 
n 

a) 

« 

PI 

rtf 

3 

0 

0 

0 

c 

1 

Ul 

1 

© 

1 

© 

-H 

-0 

<• 

ro 

4J 

© 

<  © 

0 

O' 

© 

0 

1 — 1 

c 

O  . 

0 

O 

< 

< 

a) 

0 

0* 

0 

CD 

<u 

1 

1 

Ul 

cu 

• 

< 

X 

*— 1 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
<0000000000000000000000000000000000000 
1-0000000000000000000000000000000000000 
Ul . 

Oonono«onono«o«onoftono^0fto«onor)0fi0fio^o 

•4Hri(Hnnt^rtlft-0'ONMDID(M>-OOHHM«nfJ»Vlftrt<0-ONND 


^•ooooohCD-orocMinHMDnH  _  _______ 

o«ooHxo-Hn#NOtOHHj.nort9‘**-ooMunHHn<Hnono 

1j|HONt>t®Hf)^n^nntriHOHN<n(Min<nHiDtO‘'»NOHri 

rt  n  n  rx  rx  ooO'©N-o*><orx-«0"-*fxn<*>«oN©choo~  ex  rx  n  n  n 


4ioiMf)nHNoo*<o*<noNr«n«&0'«oNa‘0*ofi-0’0 

Hn®NOOwHo.noft9-H-fflONOnHHDtfHg)ono 


~.-.-.0000000000000000000< 
I  I  I  I  I  I  I  I  I  I  I  I 


o^na)NHp)n(MD4)a)MhO'Oif)'0<0'0«>ooo‘MD-oa)(MonHNCDnHO 

oiMh-«HooMDn^'in-o-nonor>orj-i^in^HnmNO(D-«HO‘fto 

8-»Ci<r>r>^oo<n^o<^<NO-^rtf^o^<0'<©-*n<r>^)nr>Rrftf^o 
-*rtf)<o-oNOD‘OO^Hjiirtnnnnnrin(**-oo(MDMin»nw*«o 

6o66666o66^^^“«-«*«^^^«^-h-«h-«h»«666666  6  666 


■)-Mntf«-oMDtf‘0-wot«'OMDff‘OHrip)vn-oNOfrO-Nn<rn'ON 

H-H--HHHHHfii(nNrjrjrKiiMn(Knnnnnnnn 


o  0 

*-•  1 

K  111 

<  N 

cr  © 

O  ZD 

UJ  u  . 

h- 

z 

M 

n 

ui  o 

a:  1 

D  —  UJ 

©  x  h» 

©  ~  ftf 

Ul  <01 

cr  o  . 

cl  < 

I 

O' 

©  o 

H  I 

z  —  © 

Ul  X  O' 

•-*  —  h. 

o  >■  h* 

w  O  . 

II  O' 

tL  I 

III 

o 

o  o 

J-  ~  Ul 

z  x  o 

£  2? 

o  o  . 

E  -* 


-133- 


6TRATFQRD  SEPARATION  CRITERION  (LAMINAR)  F(S)  -0.02232 


o: 

>0 

>0 

\ 

o 

o 

o 

-o 

-0 

cr 

o 

o 

- 

-• 

X 

X 

Q 

o  «n 

a 

O  «n 

n 

n 

N 

O 

\ 

o  r- 

m 

m 

■* 

o  o 

- 

O  o 

u 

o 

O 

CL 

Q- 

N 

ft 

ft 

U 

c  6 

U 

6  6 

Q 

Q 

■*. 

o 

'  o  <r 

lO 

tD 

CL 

o 

«r 

0.  O  *T 

n 

n 

O 

o 

o 

u  o  o 

o 

> 

rv 

6 

o 

d  d 

o 

o 

o 

1 

xr 

«r>  o 

-0 

»n  o  «© 

> 

xr 

xr 

£L 

—  ft 

-0 

a.  ft  ft  -a 

O 

n  ft 

o  —  r*  — 

h- 

n 

n 

> 

n 

n 

a 

6  o 
i 

o 

i 

odd 

i  i  i 

d 

o 

X 

o 

xr  O 

ft 

UtOfJ 

< 

o 

o 

a: 

O'  to 

<*■ 

cr  O'  r>  «r 

o 

b~ 

'O 

Xf 

< 

o  -o 

xr 

<  o  -o  •* 

o 

U 

ft 

> 

d  ft 

|N 

a 

n 

-0 

Q 

ri 

d  ni  ri 

O 

*T 

o' 

o 

i— { 

Q 

X* 

_4 

ft 

2 

O 

■ 

a) 

< 

o  o 

>0 

<  o  o  -r 

P) 

n 

Cr> 

03 

b- 

o  o 

in 

1-0  0  ft 

X 

IO 

o 

(d 

►— 

Ui 

o  o 

pj 

lil  O  O  4 

N 

o 

o 

2 

a 

C3  -  -  . 

b- 

ft 

ft 

A. 

« 

d  d 

a-* 

O  «  o 

< 

O 

— * 

•T 

go  xr  — 

w 

w* 

a. 

— • 

-x  ft  ft 

> 

2 

b- 

C 

a 

O 

*-4 

-•  03 

'O 

^  CD  O 

m 

r> 

h- 

oi  n 

n 

ft  n  n 

b- 

n 

n 

< 

r* 

n  *o 

m  n  f> 

c r 

> 

a : 

o 

N 

< 

— *  d 

i 

0 

1 

> 

a. 

i 

o 

o 

m 

Q 

i 

1 0 

UJ 

X 

a 

0) 

2 

O  N 

ft 

OM1 

a 

— • 

< 

O  0* 

IN 

O  0*  N 

u. 

* 

o 

o 

> 

O  — 

>•  O  * 

o 

CD 

CD 

2 

E 

o 

O 

a 

6  — 

o 

o  d 

01 

M 

i  i 

X 

o 

d 

h* 

V- 

D 

o 

a 

2 

M 

UJ 

tr 

CL 

> 

'O 

>0 

H 

H  UI 

J-  UI 

H 

«-« 

— • 

(O 

ol  a 

cl  a: 

01 

E 

CD 

CD 

►< 

3 

D 

< 

O 

O 

Q 

Z  0) 

2 

Z  01  2 

a 

O 

o  to 

O 

O  01  O 

z 

o 

o 

UJ 

M  UI 

h|||h 

< 

i 

a: 

h-  <K 

»- 

K-  OC  H 

D 

<  CL 

< 

<  a.  < 

ui 

01 

Z 

<r 

z  a: 

z 

03 

o  . 

< 

o  .< 

o 

> 

— « 

w 

UJ 

<  z 

ft 

<20. 

M 

z 

IK 

£  *"• 

UJ 

H  «  UI 

H 

a. 

83  E 

CO 

CO  E  CO 

M 

0) 

u. 

O 

o 

ui 

ui 

a 

ui 

ui 

a 

a 

o 

a 

> 

-j 

*“• 

K 

03 

to 

< 

01 

CD 

i 

> 

> 

i- 

>> 

> 

r 

■f 

i 

M 

♦ 

i 

o 

z 

V 1 

-134- 
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Figure  11.-  Continued. 
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Figure  11.-  Continued 
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Figure  11.-  Continued. 
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Figure  11.-  Continued 
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Figure  11.-  Continued. 
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Figure  11.-  Concluded. 
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Figure  15.-  Measured  and  predicted  normal-force 
coefficient  on  an  ogive-cylinder 
model . 


-147- 


iment 


0.0  2.0  4.0  6.0  8.0  10.0 


x/D 


Figure  20,-  Measured  and  predicted  pressure 

distribution  on  the  windward  meridian 
of  several  ogive-cylinder  models  at 
a  =  20°  . 
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Figure  22.-  Measured  and  predicted  pressure 

distribution  on  an  ogive-cylinder 
model  at  a  =  45° . 
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Figure  23.—  Predicted  asymmetric  wake 
at  x/D  -  4.7  on  an  ogive- 
cylinder  model  at  a  =  45° 


Figure  24.-  Measured  and  predicted  velocity  field  at  x/D 
on  an  ogive-cylinder  model  at  a  =  45°. 
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Figure  25  Measured  and  predicted  velocity  field 

at  x/D  =  4.9  on  the  lee  side  of  an 
ogive-cylinder  at  a  =  22.4° 
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Figure  26.-  Measured  and  predicted  velocity 
field  at  x/D  =  4.9  on  the  lee 
side  of  an  ogive-cylinder  at 
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Figure  27.-  Measured  and  predicted  normal- 
force  coefficients  on  a  series 
of  elliptic  cross  section  bodies 


Figure  28.-  Predicted  vortex  cloud  wakes  on  a 

series  of  elliptic  cross  section  bodies 
at  x/D  =  10,  a  =  20°. 
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